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Abstract 

We build a model for the hepatic fatty acid metabolism and its metabolic and genetic regulations. 
The model has two functioning modes: synthesis and oxidation of fatty acids. We provide a sufficient 
condition (the strong lipolytic condition) for the uniqueness of its equilibrium. Under this condition, 
modifications of the glucose input produce equilibrium shifts, which are gradual changes from one 
functioning mode to the other. We also discuss the concentration variations of various metabolites 
during equilibrium shifts. The model can explain a certain amount of experimental observations, 
assess the role of poly-unsaturated fatty acids in genetic regulation, and predict the behavior of 
mutants. The analysis of the model is based on block elimination of variables and uses a modular 
decomposition of the system dictated by mathematical global univalence conditions. 

Introduction 

Metabolic and genetic model of fatty acid metabolism. Recent advances in genetics and in 
physiology show the necessity for melting together several types of cultures in order to understand 
animal and human nutrition and solve important health and economic problems. Metabolic analysis 
approaches study dynamics of biochemical pathways and employs detailed knowledge of biochemical 
reactions mechanisms [Men97, Fel97, HS96, CB95, MLR90, PSP+04, PRP+03]. Genetic functional 
studies mostly concern gene networks, recently integrating some metabolites [dJGH+04, CCRFS05, 
LSS+05]. Gene network dynamics is modeled by various methods: systems of differential equations 
[TCN03], boolean or multivalued logical automata [Kau93, STOl], Petri nets [CRRT04, MDNMOO]. 
Although virtual cell models are planned by many, present studies deal with simple cell functions. At 
a higher level of complexity, integrative approaches were applied to modeling various organs, the heart 
being one of the best studied [Nob02] . The goal of these studies is to explain physiology from molecular 
basis. The two main obstacles against this goal are the sparseness of the biological knowledge and the 
mathematical difficulty of analyzing large complex systems. 

Main pathways of carbohydrate metabolism are considered to be the cornerstones of metabolic mod- 
eling [ECB98, Cha02, TPaOO]. However, models of metabolism in eukaryotes are scarce and dedicated 
to specific metabolic pathways and organs [LK02]. Furthermore, although recent experiments pointed 
out the genetical changes induced by diet in various organisms [CWM05, ACW+03, La04, BLC+04], 
metabolic dynamical modeling is rarely considering the genetic context. Concerning mathematical com- 
plexity, one can reduce the number of numerical parameters in the models by building so-called " minimal 
models" [BIBC79, VAS+02, CWS04, TCB+06]. Even for such minimal models inferring parameters from 
data posses non-trivial problems [MSBC06] . 

In this paper we present a mixed metabolic and genetic model of regulated fatty acids metabolism 
in liver, representing a reasonable compromise among these various cultures. Parallel work by [Bel05] 
on liver focus on transport processes and do not discuss genetic regulation. The importance of genetic 
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regulation in fatty acid metabolism was recently emphasized [CF03, PLM03, DF04, Jum04]. Here; we 
trade oflF the complexity of spatial effects against the complexity of regulations and the possibility of 
analytical reasonings. Our analysis of the model do not use numerical parameters. It puts forward 
a mathematical qualitative approach that can be used more generally for the analysis of equilibria of 

complex biological systems. 

Questions raised by mixed models. Melting cultures generates many, conceptually diverse, ques- 
tions. Let us point out three such questions related to multistability, timescales and the role of genetic 
and metabolic regulations. 

Mixed gcnetical and metabolic systems are often multistable. As a well known illustration, the 
functioning of the E.Coli lactose operon is based on bistability [PJM59, YM03]. Thus, the change 
in food (lactose) induces an equilibrium switch which represents a jump from one attractor to another. 
Equilibrium switch is efficient in saving resources (enzymes) , becaiise these are produced only on demand. 
As a counterpart it is less flexible and tuning is not possible (the response is of the binary type). 

The alternative way to adapt to external changes is via equilibrium shifts. Then, during an equilibrium 
shift, an equilibrium uniqueness condition is fulfilled and there are no jumps between attractors. The 
jumps are replaced by smooth, gradual changes. 

Fatty acid metabolism in hepatocytes has two antagonistic functioning modes which could suggest 
bistability: synthesis produces fat reserves; lipolysis burns fats and produces energy. This motivates the 
first question we wish to answer: in higher organisms, does the whole of regulations produce multista- 
tionarity or a unique equilibrium of fatty acid metabolism? We shall argue that during a change in food, 
a unique equilibrium shifts smoothly between the two functioning modes and that there is no bistability. 

The second question we wish to answer is about timescales. Genes coding for enzymes need relatively 
long times to change expression levels and enzymes concentrations. On short time scales enzyme concen- 
trations can be considered to be constants. This suggests that changes of nutritional conditions induces 
processes with various timescales. We want to know whether there are any physiological consequences 
of the multiple timescales, and to find simple ways to take this into account in our modeling. For in- 
stance, fasting demands a shift from synthesis to lipolysis functioning modes. This can be done rapidly 
by metabolic control. Genetic regulation brings slower changes that push the shift further. Other slow 
processes (for instance diffusion-controlled lipid transport within different organs) could be responsible 
for other long time scales. We shall limit our analysis to only two timescales: a metabolic, fast one, and 
a genetic, slow one. 

Among the fatty acids, many work describe the interference between genes and a special class of fatty 
acids (polyunsaturated fatty acids denoted by PUFA) [Jum04]. PUFA are synthetized from essential 
fatty acids, that are taken from the diet. To the contrary, saturated and mono-unsaturated fatty acids 
(denoted by S/MU-FA) are synthesized de novo in liver. PUFA control their own oxidation as well as 
the synthesis and oxidation of the other fatty acids. The control is due to formation of complexes that 
activate or inhibit the active forms of nuclear receptors regulating the transcription of genes coding for 
enzymes involved in the corresponding pathways. 

The third question we wish to answer is about the effects of the PUFA interaction with genes, within 
normal and mutant genotypes. We shall argue that, when genetic regulations are absent (for instance, 
in PPAR- knocked out cells), the increase of PUFA concentration during fasting is stronger than in wild 
type cells. 

Main points of the paper. The paper is structured as follows. 

• Definition of several classes of partial equilibria associated with a mixed model. In 
Section 1, we introduce mixed models, that is, a differential system of equations where genetic/slow 

variables and metabolic/fast variables are distinguished. Associated to this model, we derive two 
types of partial equilibria: first, the well-known quasi- stationnarity or non- genetically regulated 
equilibrium, where genetic variables are supposed to be constant. Second, genetic partial equilibrium 
where genetic variables are supposed to be at equilibrium. Partial equilibria result from a reduction 
method consisting in successive elimination of variables of the model. 

• Construction of a mixed differential model. In Section 2 we build a differential model 
for the regulated fatty acid metabolism in liver. In order to reduce complexity, we deliberately 
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choose a simplified description of the main mctabohc pathways and of the genetic regulations. 
Although caricatural, the models exhibits the regulatory function of PUFA and agrees with recent 
experiments. Our model is quite general since we do not use explicit forms of the functions relating 
metabolic fluxes to genetic or metabolic variables. We just take into account the signs of the 
variations of elementary fluxes with respect to the variables. Equilibrium equations allow to extract 
implicit relations between variables. This approach is close to metabolic control [CB95, Fel97, HS96] 
and also to classical equilibrium thermodynamics [Cal05]. 

• Study of equilibria. In Section 3 the steady states of the differential model are studied. We find 
a sufficient condition for the uniqueness of equilibrium. This condition is expressed mathematically 
as an inequality involving partial derivatives of the fluxes with respect to metabolites (elasticities). 
We discuss a biological interpretation of this condition. 

• Qualitative validation, prediction and illustration. In Section 4 we develop the qualitative 

analysis of the model. The behavior of the model is coherent with known experimental data. More- 
over, the model has several predictions concerning the effect of suppressing the genetic regulation 
of the important nuclear factor PPAR-a, also on the role of genetic regulation for energy recovering 
at fasting. 

Using standard regulation functions, we propose an explicit set of differential equations which 
represents the dynamics of hepatic fatty acid metabolism and its regulations. We give numerical 
simulations which illustrate the main results of the paper. 

• Discussion. Section 5 is devoted to a discussion of the results and of possible extensions. 

• Mathematical method. Section 6 provides details about the mathematical method. This method 
uses a decomposition of the system into boxes or modules that are chosen in order to fulfil the 
conditions of a global univalence theorem (Gale-Nikaido). 

1 Mathematical framework: mixed model; associated partial 
equilibria states 

Mixed metabolic-genetic, slow-fast decomposition. Main carbohydrate metabolic pathways are 
rather well documented [Sal99]. Nevertheless, actual models of these metabolic pathways do not take into 

account genetic regulation. In these models [Cha02] enzyme concentrations arc parameters, rather than 
dynamical variables. We introduce here a mixed metabolic/genetic system that contains metabolites, 
proteins (products of genes), especially transcription factors and enzymes, as dynamical variables. 

Our mixed model for genetically regulated metabolism is represented by a system of differential 
equations: 



where p £ A, A being a compact subset of R*, stands for a set of external parameters. The dynamical 
variables are partitioned into two groups. Concentrations of metabolites involved in biochemical reactions 
are metabolic variables, represented by the vector X G R" . Concentrations of proteins (basically enzymes 
and transcription factors) arc genetic variables, represented by the vector Y G R™. 

Most of the genetic variables vary on timescales generally much longer than any of the metabolic 
variables: genetic variables including concentration of products of genes (enzymes, transcription factors) 
have significant variations on long (genetic) time scales tq. On short, metabolic, time scales tm << tg 
these variables can be considered to be fixed. Thus, our partition of the variables corresponds to the 
well known slow- fast decomposition of dynamical systems [Mur03]. 

In Eq.(l.l), $ is the time derivative of fast (metabolic) variables and 'S' is the time derivative of 
slow (genetic) variables. Beyond timescales, there is a physical difference between these functions. # 
is a combination of generally conservative, metabolic fluxes. * have no reason to be conservative. 
The consequences of this difference will show up in the construction of the model (in particular when 
identifying the relations among fluxes). 




(1.1) 
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Equilibrium and quasi-stationary states. Equilibria arc defined mathematically as fixed points of 
a system of differential equations and biologically as stationary states in which measurable macroscopic 
quantities stop changing. 

Given the parameter p, recall that an equilibrium state of the system (1.1) is defined by the following 
equations: 

f =*(X,Y,p)=0. ^^-^^ 

Given a decomposition of variables and a parameter p, we can define another notion of equilibrium, 
called quasi-stationary state, which is the equilibrium state of the subsystem associated to the variable 
X, constrained by fixing Y. A quasi-stationary state satisfies: 

*(X,Y,p)=0 

Y = Yo(= const.). ^ ' 

As detailed before, mixed metabolic-genetic systems are slow-fast systems. On metabolic, time scales 

tm « tg genetic variables can be; c;onsidcrcd to be fixed. Therefore, within metabolic timescale tm the 
system reaches only quasi-stationarity. After a longer time of the order tq, it reaches equilibrium. 

Slow-fast dynamics; genetically non-regulated trajectory. The slow manifold of a slow-fast 
system is defined by the equations *(X, Y,p) = 0. 

Although the definition of a quasi-stationary state (Eq.(1.3)) can always be used, this state has a 
dynamical meaning only within some conditions which are those required for the applicability of the 
Tikhonov theorem [TVS80, Was65] or of the geometrical theory of Fenichel of singular perturbations 
[Fen79]. Indeed, the slow manifold should be hyperbolically stable with respect to the constrained 
dynamics ^ = $(X, Yo,p) (see Section 5 for details). 

Under this condition, trajectories of slow-fast systems are made of two parts: a rapid part finishing 
close to the slow manifold and along which Y is practically constant and a slow part practically included 
in the slow manifold [Was65, Fen79, Mur03]. 

Let us call genetically non-regulated trajectory the set of points approximating the rapid part of 
the slow-fast trajectory, that is, following the dynamics ^ = $(X, Yo,p) and Y ~ Yq. The quasi- 
stationary state is the intersection of the genetically non-regulated trajectory and the slow manifold (see 
figure 1.1). It represents the natural intermediate stage on the way towards equilibrium. 

Genetically non-regulated model. By definition, metabolic variables on the genetically non-regulated 
trajectory are governed by a differential dynamical system for metabolic variables that we call genetically 
non-regulated model: 

rfX 

~dt 

where #g„r-(X,p) = *(X, Yo,p) are called reduced fluxes. 

Moreover, on the genetically non-regulated trajectory, the quasi-stationary equation (1.3) reduces 
to the following equations for the metabolic variables, that are called genetically non-regulated state 
equations: 

*3„,(X,p)=0 (1.5) 

Let us note that the genetically non-regulated model is the playground for classical metabolic analysis. 
Indeed, the vast majority of metabolic analysis studies do not take into account genetic regulations and 
consider that the concentrations of enzymes are fixed. 



*5„r(X,p) (1.4) 



Genetic partial equilibria. According to the classical singular perturbation theory, the constraint 
$(X, Y,p) = allows expressing the fast (metabolic) variables X as functions of the slow (genetic) 
variables Y. In this paper we perform something different. We have already introduced the quasi- 
stationarity states and the equilibrium states. Let us introduce another type of intermediate states 
obtained by equilibrating only the genetic variables. Let us call these states genetic partial equilibria. 
Contrary to quasi-stationay states (that are natural intermediate stages on the way towards equilibrium), 
genetic partial equilibria arc not dynamically reachable (because genetic variables equilibrate after and 
not before the metabolic ones). However, they represent mathematical constructions useful for the study 
of equilibria. 
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the mixed model. 


genetic partial 
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^ =*peq(X,p) 

= *(X,Ypeq(X,p),: 

*(X,Ypeq(X,p),p) =0 


^^*peq(X,p)=0 


No dynamic approxima- 
tion of the full system. 
Equilibria correspond to 
equilibrium states of the 
mixed model. 



Figure 1.1: Geometry of the mixed decomposition. For illustration we considered one metabolic variable 
X = (X) and two genetic variables Y = {Yi,Y2). The slow manifold and the genetic nuU-cline have 
the equations $(X,Y, p) = 0, and ^(X, Y,p) = 0, respectively. The slow manifold is considered 
hyperbolically stable with respect to the constrained dynamics ^ — #(X, Yo,p). Trajectories starting 
from the initial state / have a rapid part on which genetic variables are constant Y = Yq and a slow part 
in the slow manifold. The quasi-stationary state Qs{p) is the intersection of the slow manifold with the 
line Y = Yq of genetically non-regulated states. The equilibrium state Eq{p) is the intersection of the 
slow manifold with the genetic nuU-cline Y = Ypeg(X, p) which is the line of genetic partial equilibrium 
states. 



More precisely, we intend to express genetic variables as functions of the metabolic variables thanks 
to the genetic partial equilibrium equation '4'(X, Y,p) = 0. The following condition ensures that these 
equations have solutions in variables Y. 

Condition 1 (Unique genetic partial equilibrium) 

The genetic partial equilibrium equations '4'(X, Y, p) = have a unique solution in Y which is a smooth 
function of X for all X G K",p £ A. In other words, genetic nuU-clines are smooth and have unique 
intersections with the hyperplanes X — const.. 

If Condition 1 is satisfied, we denote the corresponding implicit functions Y = Ypoq(X, p). This equa- 
tion defines the genetic null-cline (see Fig. 1.1), that is, the trajectory of the dynamics ^ ~ $pcq(X,p), 
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Y = Ypcq(X, p) where $poq(X. p) = $(X. Yppq(X, p), p) are called reduced fluxes. 

On the genetic nuU-cline, we can simply express the equilibrium equations (1.2) as a set of constraints 
for the metabolic variables only. In analogy to classical thermodynamics [Cal05] we call these reduced 
constraints genetic partial equilibrium state equations for metabolic variables: 

*peq(X,p)=0 (1.6) 

The vector field #peq(X, p) defines itself a differential dynamical system governing the metabolic 
variables on the genetic nuUcline. We call this system the genetic partial equilibrium model: 

— =*peq(X,p). (1.7) 

Equilibria and equilibria shifts. Our purpose is to discuss equilibria and equilibria shifts of the full 
system. For this purpose, the partial equilibrium model is fully suitable, as follows from: 

Proposition 1.1 

The equilibria of the partial equilibrium model (Eq.{1.7)), i.e. the solutions of the genetic partial equi- 
librium state equations (1.6) are the equilibria of the full system (1.1). 

The equilibria of the genetically non-regulated model (Eq.{lA)), i.e. the solutions of the genetically 
non-regulated state equations (1.5) are the quasi-stationary states of the full system (1.1). 

State equations help us to assess the uniqueness or the multiplicity of equilibria, or to estimate the 
variations of the metabolites when the external parameters change. Furthermore, they allow to extend 
notions from the control theory of metabolism such as elasticities [CB95] to the case when genetic 
regulation is present. The essential of metabolic control is expressed by the following formula which is a 
consequence of the implicit function theorem: 

dX 
dp 

Eq.(1.8) says that the response of a metabolic system to changes of the parameters can be calculated 
from the state equation. The choice for the function ^red should be either ^gnr or ^peq depending 
on the timescale of the changes. If changes are monitored on short, metabolic timescales then one 
should consider genetically non-regulated state equations. If the changes are monitored on long, genetic 
timescales then one should consider partial equilibrium state equations. In all cases, it is important to 
compute the derivative matrix This matrix describes the resistance of the metabolic variables 

to forcings and is analogous to the matrix of elastic constants in elasticity theory [CB95]. '^^j^'' gives 
the instantaneous elastic constants, while gives the static elastic constants (resistance to slow, 

adiabatic changes, taking into account genetic readjustment). 

Comments. 

• If Condition 1 is not satisfied then there are several equilibrium state equations in the metabolic 
state variables. We shall not discuss this situation in this paper. 

• This particular type of reduction implies a certain emphasis on metabolic variables. It is justi- 
fied when genetic variables are not measured or when we are mainly interested in variations of 
metabolites. 

• The genetic partial equilibrium model represents a bad approximation for the dynamics of the full 
system. On the contrary, the genetically non-regulated model represents a good approximation 
for the rapid part of the trajectories of the full dynamical system Eq.(l.l) (see Fig. 1.1). The 
correct approximation of the slow parts of the trajectories is given by the slow/fast decomposi- 
tion, more precisely by the reduced system ^ = 'S'(Xred(Y, p), Y, p), X = Xr-ed(Y,p) where 
*(X,ed(Y,p),Y,p)=0. 



d^r 



dX 



d^ 



red 



dp 



(1.8) 
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2 Mixed model for genetically regulated fatty acid metabolism 



We provide here the details of our model and we derive the associated genetically non regulated and 
genetic partial equilibrium models. 

In different species such as chicken, rodents and humans, hepatocyte (liver) cells have the specificity 
to ensure both lipogenesis and /3-oxidation. To set ideas, all the variables of the model pertain to an 
"abstract" hepatocyte, capable of the two different functioning modes. 

2.1 Variables and fluxes for the mixed model of regulated fatty acid metabolism 
Metabolic variables. 

We have selected the most important metabolites implied in the fatty acid metabolism in liver as 
follows. Corresponding symbols for these variables are given in Table 2.1. 

• Acetyl-CoA generated in mitocondria is the first brick for building fatty acids. It is consumed in 
lipogenesis in hepatocytes, produced in oxidation. 

• Saturated and monounsaturated fatty acids (denoted by S/MU-FA) can be produced by 
the organisms from Acetyl-CoA. They can also enter the metabolism as part of the diet. 

• Exogenous polyunsaturated fatty acids (PUFA) arc implied in genetic regulations. They 
can be manufactured from essential fatty acids which can not be produced by animals. As a 
simplification, we write that PUFA can only enter the metabolism as part of the diet. 

• Energy (ATP) expresses the energy that the cell has at its disposal. 

Notice two fundamental points: first, we have introduced the level of energy of the cell as a variable. 
Second, we have divided fatty acids into two parts: the ones that are implied in the genetic regulation, 
that is PUFA, and the ones that are not. PUFA arc synthctized for essential fatty acids provided by the 
diet; hence we consider that this class can not be produced by the cell. Even if simplified, this distinction 
will allow us to model better the regulations of the metabolism. 

Parameter. The system is driven by the glucose concentration, representing food. Different nutritional 
states such as normal feeding or fasting are modeled by different values of this parameter. 

Primitive metabolic fluxes. Main metabolic processes are modeled here as primitive fluxes. They are 
represented as unstructured reactions, whose detailed mechanisms are not described. The corresponding 
symbols are given in Table 2.1. 

• Glycolysis (in which we include the Pyruvate dehydrogenation reaction) produces Acetyl-CoA 
from glucose. Glycolysis can be considered reversible. Nevertheless, we shall not study glucose 
dynamics (G is a constant); reversibility will be neither used nor rejected. 

• Krebs cycle produces energy for cellular needs from Acetyl-CoA. 

• Ketone bodies exit allows the cell to to transfer the energy stored in Acetyl-CoA to the outside; 

it represents an important source of survival during fasting or starving. 

• Lipogenesis transforms Acetyl-CoA first into citrate, then into saturated and monounsaturated 
fatty acids S/MU-FA. 

• A outtake flux allow S/MU-FA to exit liver and go to storing tissues (adipocytes). Conversely, 
the intake fiux is fed partially from diet, partially from lipolysed adipocytes. The intake flux is 
conventionally considered positive. 

• Similar intake/ outtake flux of PUFA allows PUFA to enter or exit the cell. Above diet and 
lipolysis, the intake flux of PUFA also includes a synthetic pathway consisting of desaturation and 
elongation of essential fatty acids. 

• oxidation burns all fatty acids in order to produce energy and to recover Acetyl-CoA. 
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• ATP consumption expresses the energy (ATP) the cell consumes for living. 

• Degradation of metabolites (Acetyl Co- A, S/MU-FA, PUFA) is used with a broad meaning 
including cell growth induced dilution, leaks or transfers to non-represented pathways, and effective 
molecular degradation. Negligible on the timescale of the metabolic processes, these processes can 
not be neglected on the genetic timescale. 

Functioning modes. There are two functioning antagonist modes of the fatty acid metabolism in liver: 
lipogenesis that produce reserves, fatty acid oxidation that burns reserves and produces energy. The 

choice of the functioning mode depends on nutrition conditions: a lack of food (i.e. a sustained low level 
of glucose) stimiilatcs lipolysis and oxidation; normal feed (normal glucose level) induces lipogenesis. 

Fatty acids and genetic control. A good part of the known regulation mechanisms implies transcrip- 
tion factors such as nuclear receptors PPARa (peroxisome proliferator activated receptor a) and LXRa 
(liver X receptor a). The latter is known to activate the transcription of SREBP-1 (Sterol response 
element binding protein 1) known to trans-activate different genes involved in fatty acids synthesis and 
desaturation. 

Concerning the metabolites, it has been established that fatty acids can up-rcgulatc or down-regulate 
the expression of different genes controlling their metabolism. The regulatory effect is mainly due to 
PUFA: the interaction of S/MU-FA with genes is supposed weak. More precisely, it has been proposed 

that PUFA regulates the activity of SREBP-1 and of several members of the steroid-thyroid superfamily 
of nuclear receptors such as PPARa and LXRa (for reviews see [CF03, PLM03, DF04, Jum04]). 

Genetic v£iriables. Since the genetic interactions between metabolites and fluxes is not direct, we 
consider the following abstractions for the genetic regulation variables whose corresponding symbols are 
given in Table 2.1: 

• the active form of the nuclear receptor PPAR, 

• the active form of the nuclear receptor LXR, representing in a very simplified way the regulation 
path LXRa-SREBP-1. 

• a representative abstract enzyme for each set of enzymes that are involved in S /MU-FA synthesis 
and oxidation, PUFA oxidation and ketone bodies exit respectively. Abstract enzymes production 
is controlled by LXR (modelhng the LXRa-SREBP-1 pathway) and PPAR. 



Variable (Concentration) 


Symbol 


d product 

dt 


Acetyl Co-A 


A 




Saturated and monounsaturated 


Fi 


$Fi 


fatty acids (S/MU-FA) 






Poly-unsaturated fatty acids 


F2 




(PUFA) 






Energy ATP 


T 


$T 


Active form of PPAR 


PP 




Active form of the regulation path 


L 


*2 


LXR-SREBP 






Enzymes of S/MU-FA synthesis 


El 


*3 


Enzymes of S /MU-FA oxidation 


E2 


*4 


Enzymes of PUFA oxidation 


E3 


*5 


Enzymes of Ketone body exit 


E4 


*6 



Primitive flux 


Symbol 


Glycolysis 


Gly 


Krebs cycle 


Krebs 


Ketone bodies exit 


Kout 


Lipogenesis 


Syn 


/3-oxidation of S/MU-FA 


Oxil 


/3-oxidation of PUFA 


Oxi2 


S/MU-FA fatty acids intake/outake 


Finl 


PUFA fatty acids intake/outake 


Fin2 


ATP consumption 


DegT 


Degradation of a metabolite V 


DegV 


(V = A, Fi, F2) 





Parameter 


Symbol 


Glucose 


G 



Table 2.1: Symbols for the variables, their production (expressed as time derivatives), parameter and 
primitive fluxes of the genetically regulated fatty acid metabolism. 
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2.2 Regulations for the mixed model of fatty acid metabolism 




Figure 2.1: A model for genetic regulations of fatty acid metabolism. Dashed arrows stand for genetic 
timescale actions from the origin on to target. Plain arrows stand for metabolic fluxes. Dot arrows stand 
for energetic regulations implying T. In this model, notice that a metabolite F2 (that is, polyunsaturated 
fatty acids PUFA) regulate the genetic regulators L (LXRa-SREBP-1 pathway)) and PP (PPAR-a). 



We have first defined a sketch of the model. We now intend to add the regulation relations between 
fluxes and variables. These are described by the sign of the variations of a flux when we increase the 
activity of the regulator. In mathematical terms, regulation is summarized by the set of signs of the 
partial derivatives of the primitive fluxes with respect to the variables. We consider two classes of 
regulation, namely metabolic and genetic regulations. 

Metabolic and energetic (ATP) regulations. Metabolic biochemistry has intrinsic regulation: 
substrates stimulate and products inhibit reactions. The latter effect is a consequence of the more 
general Le Chatelier principle that says briefly that effects turn against their causes. In control theory 
this means that feedback is negative, which is in favor of uniqueness and stability of equilibrium. Other 
negative feedback regulations are responses to the energetic balance (ATP/ADP or ATP/ AMP ratios) 
via direct biochemical regulation, or more complex signalling pathways [Car05]. An increase of ATP 
favours catabolic processes (lipolysis and oxidation) while a decrease favours synthetic pathways. 
We can classify metabolic regulations as following: 

• Substrate effect An increase of substrate increases the associated flux. This implies the following 
relations in the model: 

OSyn ^ dGly ^ n. dOxil ^ n dOxi2 ^ n. 5Kout ^ n, 5Krcbs ^ r\ ctDcgT ^ c?Syn ^ p. 

dA ^ ^' aa ^ ^' aFi ^ ^' ^ ^' oa ^ aA ^' aT ^ ^' bt ^ ^• 
Also, degradation reactions are modeled as DegV(V) — (5vV, (5v > where V denotes any variable 
A, Fi, F2. Hence > 0. 

• Passive or active transport effects Intake/outtake fluxes Finl and Fin2 are conventionally 
directed to the inside. Hence, they decrease when the internal concentrations of fatty acids increase: 
^ < and < 0. 

• Product negative feed-back Fluxes producing ATP are negatively controlled by ATP. Thus, 

3Gly ^ r> cjQxil , r. 30x11 , r. 3Krcbs ^ n 
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• Energy effect on fat intake Fat intake is needed to produce energy by oxidation. A drop in 
energy (ATP) stimulates fat intake. This means that < and < q. 

Genetic regulations. The role of PUFA in genetic control has been discussed in recent publications 
[PLM03, DF04, Jum04]. Although the precise mechanisms have not been proven yet, some well estab- 
lished facts can be used for modeling: 

• Fatty Acid synthesis down-regulation. PUFA inhibit lipogenesis via the LXRa-SREBP-1 
regulation path. Concerning the mechanisms of these interactions there are some hypothesis. PUFA 

could regulate the nuclear abundance of transcription factors such as SREBP-1 via the turnover 
of its mRNA and also via its proteolytic processing which is specific of SREBP's family [Jum04]. 
It has also been suggested that PUFA can bind to and modify the activity of nuclear receptors 
PPAR and LXR: active PPAR and LXR are heterodimers with RXR (retinoid X receptor): PUFA 
could prevent nuclear receptor LXR from forming a heterodimer with RXR, and therefore blocks 
its activity as a transcription factor [Jum04].. 

• Oxidation up-regulation. PUFA stimulate their oxidation as well as the oxidation of Fi since 
they activate PPAR. The detailed mechanism is not known: it either cooperative stimulation of 
transcriptional effect of PPAR or indirect (active PPAR is a heterodimer with RXR; preventing 
LXR/RXR formation fatty acids shifts the equilibrium toward PPAR/RXR formation). 

• Ketone exit up-regulation. The mitochondrial HMG-CoA synthase, a key enzyme of the ketone 
body formation is known to be transactivated by PPARa; in vivo PPARa activation leads to an 
increase of ketone bodies exit [La04]. 

We translate this biological information into several relations between the variables. Notice that since 
these regulations imply genetic interactions, they occur only on long (genetic) timescale tq- 

• PUFA (F2) activates PPAR (PP) and inhibits active-LXR and SREBP-1 (L): II7 > 0, ||f < 0. 

• LXR and SREBP-1 (L) triggers Ei production (where Ei models S/MU-FA synthesis enzymes): 

• PPAR (PP) triggers the production of E2 (S/MU-FA oxydation enzymes), E3 (PUFA oxydation 
enzymes) and E4 (ketone exit enzymes): > 0, |pp > 0, |pp > 0. 

• Degradation effects occurs on each genetic variable. 

• Abstract enzymes Ei stimulate the corresponding fluxes. > 0, > 0, > 0, > 
0. 

2.3 Differential model for the regulated fatty acid metabolism 

The graphical representation of the model is shown in Fig. 2.1. Table 2.2 summarizes a differential model 
including the above described relations among metabolic and genetic variables. Let us summarize the 
notations: 

• X = (A, Fi, F2, T) is the set of metabolic variables; 

• Y = (PP, L, El, E2, E3, E4) is the set of genetic variables; 

• p = G is the parameter of the model; we suppose that it takes values inside a compact interval, 
Gg[0,G 

maxj * 

• $ : ^ and * : M;^° ^ are defined such that: 

r f = $(X,Y,p) 
I f = *(X,Y,p) • 

The differential model was built as follows: 
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• ^ = Gly(G,T) + niOxil(Fi,T,E2) + n20xi2(F2,T,E3) 

-Krebs(A, T) - Kout(A, E4) - niSyn(A, T, Ei) - 5aA 
^ = Syn{A,T,Ei)-Oxil(Fi,T,E2) + Finl(Fi,T)-aF^Fi 
^ = -Oxi2{F2, T, E3) + Fm2(F2, T) - 5f,F2 
^ = QGGly(G, T) + ttKKrcbs(A, T) + QoiOxil(Fi , T. E2) 
+ao2 0xi2(F2, T, E3) - asSyn(A, T, Ei) - DegT(T) 
' ^ = *i(F2) - 5ppPP 
^ = 4'2(F2) - 5lL 
^ = *3(L)-5eiEi 

# = 5'4(PP) - 5b, E2 

# = *5(PP) - 5B3E3 

[ ^ = xi/s(pp) - 5e,E4 



tlux 
;» va r. 
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Kout 
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+ 






































G 


+ 













































Table 2.2: Differential model for the regulated fatty acid metabolism: equations and constraints. The 
flux of each metabolic variable is obtained as a mass balance of primitive fluxes. The dependance of 
primitive metabolic fluxes and of genetic production terms on variables is given in the table above. 

• The production $a, ^Fi, ^Fa; of each metabolic variable is obtained as the sum of primitive 
fluxes that produce or consume the metabolite. 

• Primitive fluxes arc treated as single reactions with simple stoechiometry. Thus, the fluxes Gly, 
Krebs, Oxil, Oxi2, Syn are considered to have the stoechiometries G ^ A + acT, A UkT, 
Fi — > mA + aoiT, F2 —>■ niA + ao2T, niA + asT — > Fi respectively. 

• Degradation reactions of metabolites are supposed to be linear: DegV(V) = SyV where V denotes 
any variable A, Fi, F2. 

• The fimctions expressing variations of the genetic variables (PP, L, Ei, E2, E3, E4) were not 
detailed because mechanisms arc still unknown. Instead, each function 'I'; has been decomposed 
into a non-negative production term i and a linear degradation term. Qualitative information on 
the production term is translated into a set of signs of its partial derivatives as in Table 2.2; for 
the time being these are the only a priori constraints to the model. 

2.4 Reduced models for the metabolic veiriables 

Let us construct the two reduced models and state equations for the metabolic variables: flrst, the genetic 
partial equilibrium model and its associated state equations; second, the genetically non-regulated model 
and its state equations. 

In order to obtain the genetic partial equilibrium model we must eliminate the genetic variables from 
their equilibrium equations, i.e. we must solve the subsystem (genetic partial equilibrium): 



Notice that in the subsystem (2.1), the global variable F2 is a parameter. If a solution of this 
subsystem exists, then the genetic variables PP, L, Ei, E2, E3, E4 are expressible as functions of F2. It 
can be easily shown that this is possible in a unique way, meaning that the Unique Genetic Partial 
Equilibrium Condition 1 is fulfllled: 




*4(PP) - <5b2E2 =0 
*5(PP) - fesEa = 
*6(PP) - <5e4E4 =0. 



(2.1) 
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Gly(G, T) + niOxilpcq,9nr(Fi, Fa, T) + n20xi2poq,g„r{F2, T) 
-Krebs(A,T) - Koutpeq,3„r(A, F2) - niSynp^^ ^^^(A, F2, T) - 5aA f[ 

Synpeq.gnr { A, Fa , T) - Oxilpeq.gnr- (Fl , F2 , T) +' Finl (Fj , T) - F 

-Oxi2peq,9„r(F2,T) +Fm2(F2,T) - Sp^Fa 
«GGly(G, T) + aKKrebs(A, T) + aoiOxilpoq,9nr(Fi, Fa, T) 
+a020xi2peq,9„,(F2, T) - asSynp,q,g„,(A, F2, T) - DegT(T) 




a) genetically non-regulated model. 



a Mux 
a variahin 



gnr 
peq 



Gly Krebs Kout Syn Oxil Oxi2 Finl Fin2 DcgT 



+ 




+ 


+ 











+ 











+ 




+ 





+ 















b) genetic partial equilibrium model. 







Table 2.3: Reduced models for the metabolic variables of the regulated fatty acid metabolism. The 
reduced primitive fluxes and their regulations are also represented graphically. 



Proposition 2.1 

For any nonnegative value of F2, the partial equilibrium equations of genetic variables (Eq.(2.1)) admit 
a unique solution, i.e. the Condition 1 is fulRlled. 

Proof. Solving the partial equilibrium equations gives readily the values of genetic variables: PPpcq(F2) = 

S^M'P2), Lpeq(Fa) = ^V^2(F2), Elp,q(F2) = 3l-^3(Lpeq(Fa)), E2peq(F2) = ^V^4(PPpeq(Fa)), E3peq(Fa) = 

^^5(PPpcq(Fa)) and E4peq(Fa) = 3^^^6(PPpcq(Fa)). ■ 
Proposition 2.2 

At genetic partial equilibrium, the derivatives of the values of the genetic variables with respect to F2 
satisfy 

'^PPpcq ^ r, I^Lpeq C^Elpcq C^Eapcq I^Espcq rfE4pcq 

dFa ' dFa ' dFa ' dFa ' dFa ' dFa 

Proof. We use the chain rule formula. For example: ''^p^"'' ~ s^Wi ^- " 

The next step of the reduction is to express the primitive fluxes at partial equilibrium as functions 
of metabolic variables only. For instance, Oxi2pcq(F2, T) = Oxi2(F2, T, E3p(,q(F2)). 

Proposition 2.3 

At genetic partial equilibrium, the fluxes Syn, Oxil, Oxi2 and Kout become functions of A,Fi,F2,T 
only, denoted by Synp^q, Oxilpeq, Oxi2peq, Koutpeq. The dependence of these functions on F2 satisfy: 

^SyHp^q aOxilpcq 90xi2pcq gKoutpeq 

aFa ' dF2 ' 9Fa ' ^Fa 



Proof. By definition Synp,q(A,T) = Syn(A, T, Eip,q(F2)). Then = By Lemma 2.2 

and Table 2.1 it follows that ^^f^ < 0. 

The sign of the other derivatives is computed in a similar way. ■ 

In genetically non- regulated states genetic variables are constant, equal to their initial values: 
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Ei(t) = Ei(0), E2(t) = E2(0), E3(t) = E3(0), E^it) = E^iO) (2.2) 

Primitive fliixes are functions of metabolic variables only and their derivatives follow directly from 
the Table 2.1: 

Proposition 2.4 

In genetically non-rcgulatcd states, the fluxes Syn, Oxil, Oxi2 and Kout become functions of A, Fi, F2, T 
only, denoted by Syn^^^, Oxil gnr, Oxi2gnr, Koutg„r- The dependence of these functions on F2 is as 
following: 

^^y^gnr r» dOxilgnr r» 90xi2gnr r\ dKoutgnr n 

Table 2.3 summarizes the reduced models in the two situations together with the table of constraints 
(signs of the partial derivatives of the reduced primitive fluxes with respect to the metabolic variables). 

3 Equilibrium and quasi-stationary states: existence and unique- 
ness 

The process of elimination of variables simplifies the study of equilibrium states of the model. We have 
eliminated the genetic variables in order to obtain the state equations for metabolic variables. These 
state equations have to be solved. The solutions of state equations are equilibria or quasi-stationary 
states. In this section we focus on the number of solutions. As discussed in the introduction if equilibria 
are unique, modifications of the variables induced by (slow) changes of the parameter of the model (food) 
are smooth equilibrium shifts. The same is true for quick changes if quasi-stationary states are unique. 
In order to find solutions of the state equations we proceed by further elimination of variables. The 
technical details are given in Section 4. 

3.1 Existence of an equilibrium and of an quasi-stationeiry state 

The concentration of metabolites results from the balance of production fluxes and degradation or con- 
sumption fluxes. 

Until now, we have assumed the following conditions on the elementary fluxes: 

Condition 2 (Flux global constraints) 

• The fluxes are differentiable functions of the concentrations and satisfy differential constraints 
(signs of partial derivatives), summarized in Table 2.2. 

• Degradation terms are linear: DegV(V) = SyV where V denotes any variable P, A, Fi, F2. 

These hypotheses are very mild. In order to go on with the analysis, we now add some more assump- 
tions which are natural and not restrictive: 

Condition 3 (Boundary and asymptotic conditions) 

• In the absence of substrates all fluxes vanish. 

• All fluxes except degradation saturate at high concentrations of metabolites. 

• ATP consumption is an increasing function of ATP with no saturation effect, that is limT_»oo DegT 
= +00. This is consistent with the fact that cells can not store ATP. 

• There exists a recovery effect on each metabolic variable. By recovery effect we mean that if 
a variable is zero, then at least one elementary flux that produces the variable is activated. In 
particular, if the cell contains no PUFA, then PUFA enter the cell. 

By a mathematical argument, we can prove that under these assumptions, our model has at least a 

quasi-stationary state and at least an equilibrium. Details arc given in Section 6. Let us suppose that 
the glucose concentration can change between and a maximal value Gmax- 
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Theorem 3.1 (Existence of equilibrium) 

Let us suppose that the conditions 2 and 3 are satisfied for all glucose concentrations within an interval 
< G < Gmax- Then, the genetically non-regulated model (Eq.{2.2) and Table 2.3) and the genetic 
partial equilibrium model (Eg. (2.1) and Table 2.3) for fatty acid metabolism admit at least an equilibrium 
state for every < G < Gmax ■ 

According to Proposition 1.1, we derive a result about the fuU model. 
Biological prediction 1 

Under the conditions of Theorem 3.1 the regulated fatty acids metabolism model described in Table 2.2 
has at least a quasi-stationary state and at least an equilibrium. 



3.2 Uniqueness of equilibrium 

Various methods using the interaction graph provide sufficient criteria for the uniqueness of equilibrium 
of a differential model [Tho81, Sno98, Gou98, Sou03]. In this section, we show that those methods do 
not apply to our case and we propose a different method. 

Interaction graphs for metabolic variables. Oriented interaction graphs can be defined for systems 
of differential equations [Sou03]. The interaction graph gathers information relative to the (direct or 
indirect) action of a variable on another one and is thus important in the theory of response [RLS"*'06]. 
The interaction graph can be computed at genetic partial equilibrium or at fixed genetic variables (in 

genetically non-regulated states) as follows. There is an arc from Xi to Xj whenever ^ meaning 

that Xi has an influence on the flux of Xj . The sign of the regulation arc is the sign of the derivative . 
At fixed genetic variables the interaction graph gathers purely metabolic infiuences between metabolites; 
it corresponds to response on timescales that are too short to allow for genetic readjustments. At genetic 
partial equilibrium the interaction graph gathers both metabolic and genetically mediated influences 
between metabolites. 

In order to build the interaction graphs, we need to compute the signs of the derivatives of the fluxes 
with respect to the metabolites. This is done in the Table 3.1 by using Table 2.3. The sign of the fluxes 
of metabolites is computed like in the following example (corresponding to genetic partial equilibrium): 



a*? 



dFi 



J-[Synp,q(A, F2, T) - Oxilpeq(Fi, F2, T) + Finl(Fi, T) - ^f^Fi] 

- — Oxilpeq(Fi,F2,T)+ — Finl(Fi,T)- 6f, < 0. 
Or 1 Or 1 



>0 



<0 



>0 



a flux 
r? variable 






$F2 'i'T 


A 




+ 


(/) 


Fi 


+ 




+ 




+ 







p gni- 
r 2 

peq 


(11) 




+ 


T 


{III) 


(IV) 


iv) 


G 


+ 





+ 




a) non-gcnetically regulated 
model (flxed genetic 
variables, metabolic 
influences) 




' * irih 


V 













b) genetic partial 
equilibrium model (genetic 
and metabolic influences) 



Table 3.1: Signs of the partial derivatives of the fluxes and interaction graphs for the models derived 
from the regulated fatty acid metabolism model at fixed genetic variables and genetic partial equilibrium. 
The roman numerals indicate partial derivatives and regulation arcs whose signs depend on the point 
where they are computed. 
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Simple topological conditions for the uniqueness of equilibrium do not apply to any of these interaction 
graphs. Thomas rule asserts that the absence of positive loops in the interaction graph for all values 
of the node variables and external parameters is a sufficient condition for the uniqueness of equilibrium 
[Tho81, Sno98, Gou98, Sou03]. Unfortunately, it appears that Thomas rule can not be applied in our 
case as stated in the following proposition. 

Proposition 3.1 

The interaction graphs of the genetically non-regulated and of the genetic partial equilibrium models for 
the regulated fatty acid metabolism have both at least a positive loop for any value of the variables and 
of the external parameter G. 

Proof. There exists a positive loop between A and Fi (see Table 3.1). ■ 

In the following, we develop another method to prove the uniqueness of equilibrium. This method 
proceeds by successive eliminations of variables until we are left with only one equation. Then we write 
down a sufficient condition for having a unique solution to that last equation. The success of the method 
depends on the appropriate gathering of variables that can be simultaneously eliminated in an unique 
way. This is described in full generality in Section 6. 

Uniqueness of equilibrium and of the quasi-stationary state A maximal set of variables that can 
be eliminated in a unique way from the state equations is {A,Fi,F2}. In order to eliminate them, we 
consider the following state equations, at fixed T and G (those are both considered as parameters here): 

$^^<^(G,A,Fi,F2,T) = 0, (A,Fi,F2,T) = 0, (F2,T) = 0. (3.1) 

Let us recall that is cither ^gnr or $peq depending on the situation: at fixed genetic variables 
or at genetic partial equilibrium. 

Proposition 3.2 

Suppose that the Conditions 2 and 3 are fulhlled for < G < G^ax and T > 0. TJien the system (3.1) 
admits a unique solution for any given pair of values (G,T) G [O,^^^^] x l^+> both at fixed genetic 

variables and at genetic partial equilibrium. 

The functions A^^\G, T), Fi(^)(G, T), F2(^)(G, T) expressing this solution are diffcrentiable in (G, T). 

To ensure uniqueness of equilibrium, we need to eliminate T from the state equation 
$W(rp,G) = $5f'^(G, A(i)(G,T),Fi(i)(G,T),F2(^Hg,T),T) = 0. 

A sufficient condition for the unique elimination of T is < 0. In Section 6 we show that —^f- is 

a linear combination with positive coefficients of derivatives with respect to T; this allows us to restate 
the uniqueness conditions as follows: 

Condition 4 (Strong lipolytic condition) 

The following inequality is satisfied for any (G, T) e [0, Gmax] x 

aA(^> aKrebs ok 9Syn-, , aFi'^) „, aOxil , dFj'^f^, dOxil aSyn , ^, 90302 \ , d^T ^ n 

-gT^(Q!K-gA ^"Sa ) + -aT^^Oi-gpT + -OT^l^Oi-gpT ~ "S^-gFT + "02-55;-) + -g^ < 0. 

Theorem 3.2 (Uniqueness of equilibrium) 

Suppose that the Conditions 2, 3, 4 are fulfilled at fixed genetic variables and at genetic partial equilib- 
rium, for every G E [0, Gmax]- 

Then, the state equations (3.1) for metabolic variables admit a unique solution for every G G 
[0, Gmax] , both at fixed genetic variables and at genetic partial equilibrium. The concentration of metabo- 
lites and ATP at equilibrium are diffcrentiable functions ofG: A^^) (G),Ff ^ (G),F^^^ (G),T(2) (G). 

In other words, the equilibrium and the quasi-stationary states of the full model for genetic regulations 
of fatty acid metabolism are unique. 

Biological prediction 2 (Shift mechanism) 

When the Conditions 1,2,3,4 are satisfied, the mechanism allowing the change of functioning modes 
(lipogenesis and oxidation) is a shift, both for rapid response (quasi-stationarity) and for slow static 
response (equilibrium). 
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3.3 Computable version and biological significance of the strong lipolytic 
condition 

In practice, the strong lipolytic condition 4 is not rc^adily exploitable because it contains the derivatives 
of implicit functions A^^^, Fi^^-*, F2^^-'. In order to give a computable version of it, we express the above 
derivatives by using control coefficients and metabolic elasticities. 

Control coefficients, Metabolic elasticities. In metabolic control language, a control coefficient 
quantifies the dependency of a fltix on an enzyme activity. The general idea of control coefficients is 
that they permit to compare the strength of fluxes variations one with respect to the other. Following 
[CB95], we call control coefficient of a flux the derivative of the logarithm of the flux with respect to 
the logarithm of the enzyme concentration. In our problem, although they are not enzymes, F2 and T 
play regulatory roles on the fluxes through genetic and metabolic factors. Notice that these effects are 
absent in classical metabolic control analysis (this focuses on the effect of enzymes). Thus we choose to 
call non-logarithmic control coefficients all the following quantities: 

• genetic (non-logarithmic) control coefficients: 



DOxi2 
itp. 



aoxi2B 



dF2 



nOxil 
Up 



dF2 



R 



Syn 



d¥2 



9Koutp 



F2 — 



dF2 



ATP (non- logarithmic) control coefficients: 



oFini _ SFinl 



DOxi2 



DFin2 
-TLX 



90xi2p 



dT 



' dT ' 



T-»Krebs 



dSyn 



dT 
9Krebs 



dT 



oOxil 
-fix 



aoxilp 



dT 



pGly 



dGly 
dT ' 



These quantities were defined such that they are all positive (see Proposition 3.3). This sign choice 
simplifies the identification of balances in the interaction graph. For instance the interactions in the 
interaction graphs with undetermined signs can be expressed as sums and differences of positive control 
coefficients: 



(///) 



oKout 



= R!^ 



{IV) 
(V) 



R, 



Syn _|_ ^Oxil 



oOxi2 

-Tlrp 



R^ 



R'. 



The non-logarithmic control coefficients can be defined for all values of the metabolic variables. 
In particular, at fixed genetic variables the relations {R'p^^^)gnr = 0, {F^^)gnr = 0> (-RF2°"*)grar = 
express the absence of genetic regulation. The coefficient J?p^'^ is the sum of a positive metabolic 
and a positive genetic contribution. The genetic contribution vanishes at fixed genetic variables, hence 

gnr < (-RF^'^)peq- 

Following [CB95], we call elasticity the derivative of the logarithm of the rate with respect to the 

logarithm of the substrate concentration. In our setting we call non-log a,rithmic elasticities the following 
coefficients: they quantify how rates and fluxes of a metabolite depend on this metabolite. Recall that 
the rates $y are sums of primitive fluxes in Table 2.3. 



tot 5$Fi 

Xfi = -- 



9Fi 



tot ^ _^5f2 



tot _ d^A 

XA qj^ , 



dA ' 



Krebs 

Xa 



SKrebs 
dA ' 



Oxii 90xilp 
Xfi = 



dFi 



Then the last interaction with undetermined sign in the interaction graph can be expressed as (/) : 

d^T _ „, ...Krcbs 

-or - "KXa 



Syn 
"SXa • 



It will be useful in the following to introduce the following elasticity ratios which evaluate the con- 
tribution of one primitive flux to the total elasticity of a metabolite. 
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, Oxil Syn 
Oxil _ XFi Syn _ ^IXa 

^Fi — ^tot ' Pa - tot ■ 
Xfi Xa 

Also it will be useful to consider the genetic control coefficient -Rp^'^ as an elasticity and to define 
its associated elasticity ratio: 

7jOxi2 
Oxi2 _ ^F2 
^F2 - ^tot ■ 
XF2 

The values of the above defined quantities must satisfy the following: 
Proposition 3.3 

The genetic control cocfEcients, the ATP control cocfEcicnts and the elasticities are non negative. The 
elasticity ratios are well defined, non negative and less than 1. 

Proof. From Table 2.3, we can easily prove that genetic control coefficients and ATP control coefficients 
are all positive. At fixed genetic variables , some of the coefficients vanish (i?F^'^ = Rp^^'^ = i?F°"* = 0). 
Prom this analysis and Table 3.1, we get that all elasticities and the implicit ATP control coefficient are 
positive and the elasticity ratio are well defined, non negative and smaller than 1. ■ 

Let us now restate the results of the previous sub-section in terms of control coefficients and elastic- 
ities. First, let us suppose that the stocchiometry condition is fulfilled: 

Condition 5 (Stoechiometric condition) 

The stoechiometric coefficients satisfy the inequalities: as < ctoi < niac, ?^2Q!ol < niao2- 

The stoechiometric condition can be checked from biochemical data. Indeed rii , n2 are the numbers 

of Acctyl-coA and aoi,ao2 arc the numbers of ATP molecules produced (in the average by different 
fatty acids of the same type) by oxidation of a molecule of S/MU-FA, or PUFA respectively, as is the 
average number of ATP necessary for the synthesis of a molecule of S/MU-FA. /3-oxidation produces 5 
molecules of ATP for each released molecule of Acctyl-coA, thus aoi = 5(ni — 1), ao2 = 5(n2 — 1). 
PUFA have in the average longer chains than de novo synthesized fatty acids, meaning that n2 > ni. 
We deduce n2Q;oi < n\ao2- Synthesis consumes less ATP than oxidation (for example, for palmitic acid 
as = 23, aoi — 35); by generalization, we get as < aoi- Finally, aQ = 7, aK — 12 (these represent 
the number of ATP molecules produced by glycolysis and Krebs cycle per each molecule of Acetyl- Co A), 
from which aoi < »^iQ;g- 

Let us define the following combinations of control coefficients: 
A = XpOf'i, X = nr{ao^-aspT "rn^aKpT"!: 

B = Bii?irM°*+s2i?fr/XF2 +s3<vxr.*+i34PFf', 

= X-n,{aoi-as){l-p^rp^f), B^ = aoi{l - pf'^p^f^) - X/n,, 

= X(l-pO-i), Bi = mao2{l-p''rpdf')+n2/niX-n2aoi{l-ip''rfpdf'), 

?Oxil 
T 

[n2/niX + n2asp^/''{l - p^f^) + {um - n2aoi)(l - Px''' pdf^)]R?"^ + 
-{as + X/m)p^f' + p^r]B^^^''^ + [mias - aoi) + X]{1 - p^f')Rp^", 

D = [X/m - aoi(l - p^r)]RT'''''' + n2aspf''{l - p^f^)Bg^'^ + ni(aoi - as){l - pT)P°f'RT- 

(3.2) 

The strong lipolytic condition is equivalent to the following, more explicit condition: 
Proposition 3.4 

Tiie strong lipolytic response condition reads: 

7l(i?5'"i - R^""'^) + B(ii^™2 _ ^Oxi2^ +OD (3.3) 

Furthermore, if the stoechiometric condition 5 is fulfilled, then the combinations of control coefRcients 
X, A, Bi,Bi, C, D defined by Eq.3.2 are positive. 



C = [X/ni + {niaG-aoi){l-p^/''p^f)RT'' + [(Xoi+niaK+XR^ 
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In fact, fluxes of are much smaller than fluxes of F^. Further simplification of the condition is 

reasonable: 

Proposition 3.5 

If - iJ?'^'^)! «< - B^''^ then the strong Upolytic response condition reads: 

A{R^"'^ - ijg'^'i) + OD where A, C,D>0. (3.4) 

Comments. 

• In Eqs.(3.3),(3.4) all control coeflicients are functions of G and T because the strong lipolytic 
condition 4 has to be checked for all (G, T) G [0, Gmax] x R+- 

• The conditions in Eqs.(3.3),(3.4) may seem more complicated than the Condition 4. Nevertheless, 
they are readily computable and they provide the biological significance of the strong lipolytic 
condition. Eq.(3.4) is fulfilled if B^™^ — R^^^^ is large enough, which means that the energy 
variation has a suflaciently strong eff'ect on the arrival of fatty acids inside the cell. This justifles 
the name strong lipolytic response condition. 

Biological prediction 3 

Under the hypothesis of 3.5, the strong hpolytic condition means that the energy variation has a sufE- 
ciently strong effect on the arrival of fatty acids inside the cell. 

4 Validation, prediction and illustration of the model 

Our model of regulated fatty acid metabolism can be considered at difl'erent levels. 

At a qualitative level, the model consists of a set of differential equations, together with the differential 
constraints in Table 2.2. The functions giving the fluxes are not specified, neither the numerical constants 
involved in these functions. In order to validate the model or make some predictions, we provide sufficient 
qualitative conditions under which the model has a certain behavior. 

If the behavior has been proven experimentally, the conditions should be added to the qualitative 
model as extra constraints to provide a valid model. If the behavior is a hypothesis not yet proven, 
the satisfiability of the sufficient condition provides predictive models. For instance the strong lipolytic 
condition guarantees the uniqueness of equilibrium. We do not know experimentally whether the equi- 
librium is unique or not. This could be tested by response experiments, by looking at the absence or at 
the presence of (as in the case of the operon lactose of E.coli) hysteresis in the response curves. 

In this section, our aim is to determine suflicient conditions either to render from the known behaviors 
of the system or to have predictions on this behavior. Ideally, the conditions should accept biological 
interpretation. The type of behaviors that we intend to discuss within this approach are about signs 
of variations of metabolite concentrations in fasting/refeeding protocols (rather standard in biological 
studies of metabolism). 

Quantitative versions of our model can be obtained by replacing the undetermined functions in 
Table 2.2 by specific functions containing numerical constants. The advantage is that we can probe the 
dynamical behavior, much easier than in qualitative models. Specifying realistic functions and constants 
is an enormous task for such complex biological systems. In vitro measurements of the kinetical constants 
are rarely available and are not always reliable. Furthermore, low complexity abstractions are only very 
approximate models. Thus, only robust features of dynamics of the model (that are stable against changes 
of the parameters or of the forms of functions) are meaningful. We shall use quantitative versions of the 
model as illustrations of robust dynamical behaviors. 

4.1 Fatty acids concentration increase at fasting 

Let us suppose that fluxes in the qualitative model satisfy the Conditions 2,3,4. Then, a unique 
equilibrium state and a unique quasistationary state exist for any value of the glucose concentration 

max] ' 

The next result is about the dependance of PUFA concentration on glucose. In order to state our 
result, let us notice that at equilibrium or at quasi-stationarity the control coefficients i?^™^ , i?^'''^ 
depend only on G. 
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Proposition 4.1 

This proposition can be stated in biological terms as follows. 



Suppose that the Conditions 1-5 are satisfied. Let (G) be the value of T at equihbrium or quasi- 
stationarity (function ofG). Then the sign of is equal to the sign ofB^''''^ - ii?'"^ 



Biological prediction 4 (response of PUFAs during fasting) 

Suppose that the Conditions 1-5 are satisfied. The following predictions are valid for rapid (at quasi- 
stationarity) as well as for slow (at equilibrium) response: 

• If (i?:^:^"^ _^Oxi2-)^^^^ > for any < G < Gmax, then PUFAs increase during fasting and 
decrease during feeding. 



If {R^'''^ -^T'''^)egqs < for aiij < G < Gmax, then PUFAs decrease during fasting and 
increase during feeding. 



In this case, the vahic of i?^'"^ — i?!^'^'^ imphes two distinct behaviors. This means that the quahtative 
constraints associated to the differential are not precise enough to decide which behavior occurs. However, 
in this case, some biological information is avalaible about the behavior of PUFAs during fasting. We 
will use this information to refine the model with an additional qualitative constraints. 

More precisely, Lee ct al. [La04] studied for wild-type and PPAR-/- mutant murine liver, the fatty 
acids profiles in triglycerides (TG), which are the predominant (> 50%) hepatic fatty acids and also in 
phospholipids (PL) which go into celhdar membranes. Let us recall that TG and PL are storage forms of 
fatty acids and that PL contribute much less than TG to the total fatty acid mass. These authors [La04] 
show that for wild type hepatocytes after 72h of fasting fatty acids profiles do not change significantly in 
PL, but there is a strong increase of TG and of their fatty acids constituents, in particular PUFA. Based 
on these experimental findings, we make the hypothesis of a mass increase diiring fasting, of regulating 
PUFA in the hepatic cell, and look for sufficient conditions ensuring this behavior in our qualitative 
model. This is consistent with the regulation role of PUFA: in order to trigger oxidation a persistent 
increase of PUFA concentration is needed inside the cell. A qualitative reasoning using an extended 
genetically regulated model of lipogcnesis [RLS+06] gives further support to this hypothesis. 

Consequently, we have to add — i?ij?^'^)^^ > as a qualitative constraint to the model to fit 

with behaviors observed in [La04] . By this way, we make a refinement of the model. 

Condition 6 (Experimental constraint added to the model) 

A necessary condition for the observed behavior of PUFAs at fasting (increase) is (i?^'"^ — i?^'^'^)^^ > 
0, meaning that at equilibrium (or quasi-stationarity) the intake control overcomes the oxidation control 
for PUFA. 

Remark The strong lipolytic response condition, which is sufficient for the uniqueness of equilibrium asks 
that the intake control overcomes the oxidation control for S/MU-FA fatty acids (the major contribution 
to the total mass of fatty acids). Fasting experiments imply that the strong lipolytic is satisfied for 
all fatty acids, including PUFA. Nevertheless, this is not an experimental proof for the uniqueness of 
equilibrium, because one needs the strong lipolytic condition to be satisfied, not only at equilibrium (or 
quasi-stationarity), but in all T constrained states as well. Unfortunately, these states are not accessible 
experimentally. 

The next result is a dynamical one. This may seem paradoxical, because in this paper wc study equi- 
libria and equilibria shifts which is a statical problem. Nevertheless, the definition of quasi-stationarity is 
based on dynamical timescales. Since quasi-stationarity is reached before equilibrium, comparing these 
two states provides some information on the dynamics of the mixed differential model. 

In order to formulate the next result let us denote by i?eg, Bqs the values at equilibrium and at 
quasi-stationarity of the combination of control coefficients B by Eq.(3.2). 



Proposition 4.2 

If the Conditions 1-6 are satisfied and if furthermore Bf.g > Bqs then 



dG 



> 

qs 



dF2 



(2) 



dG 



eg 
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This result nicians that when G clccrcasos (fasting), the vahie of F2 at quasi-stationarity is greater 
than the value at equilibrium. Thus, we can state the following prediction of the model. 

Biological prediction 5 (Overshoot of fatty acid concentration) 

Under the hypothesis of Prop. 4.2, the curves representing PUFA concentration during fasting must 
show an overshoot: the increase in concentration is greater immediately at quasi-stationarity than later 
at equilibrium. 

Comments. 

• The condition Beg > S,, is equivalent to Bi {R§^'')^^+B2 {Rfr')eg+^^ (^Fr)^/^4^|ff > 

0, with Bi > 0, B4, ^Q^f' > 0. This means that even if B2, -B3 are negative the oxidation control 
term is strong enough to win. At fasting, this is a plausible supposition. 

• The existence of an overshoot of faty acids concentration during fasting depends on the dynamical 
accesibility of the quasi-stationary state. This state is accesible for a discontinous step-like glucose 
input, but may not be accesible if the glucose input drops slowly (see also the sections 4.4, 6.1). 



4.2 Genetic regulation reinforces energy homeostasis 

With the same methods, we prove in Section 6 how the values at equilibrium respond to the variations 
of the entering node G: not surprisingly, we predict that ATP decreases at fasting. 

Proposition 4.3 

If the Conditions 1-5 are satisfied, then ^^g^ > at at equilibrium or quasi-stationarity, where T(^)(G) 
is the value of T at equilibrium or quasi-stationarity (function of G). 

In biological terms, comparing the value of at quasi-stationarity (representing a model with no 
genetic regulation) and at equilibrium allows understanding the role of genetic regulations. 

Biological prediction 6 (response of ATP during fasting) 

If the Conditions 1-5 arc satisfied, then ATP decreases during fasting and increases during feeding. This 
prediction is valid for rapid (at quasi-stationarity) as well as for slow (at equilibrium) response. 

Proposition 4.4 

Under the hypothesis of Prop. 4.2, then (^^) > (^^) > 0. 

V / qs V / eq 

Let us emphasize that the derivative quantifies the energy buffering effect: the lower is this 

derivative, hence the lower is the variation of T for a fixed variation of G, the stronger is the energy 
buffering effect. Thus, we can state the following: 

Biological prediction 7 (Role of genetic regulations in energy homeostasis) 

Under the hypothesis of Prop. 4.2, genetic regulation reinforces the energy buffering effect. 



4.3 PPAR knock-out reduces energy buffering and increcises PUFA at fasting 

Our qualitative model can predict the effects of various genetic perturbations. Let us consider here the 
effect of PPAR knock-out. Without PPAR, there is no longer a genetic control on oxidation, therefore 
we expect to have less energy buffering on fasting. Less obvious is what happens to the concentration 
of PUFA. We can predict the behavior of PUFA. under the same hypothesis as above. Let F2*^^-'(G) be 
the value of PUFA concentration as a function of G. Also, let BwT,eq, BppAR-/-,eq be the values at 
equilibrium in wild type and mutants of the coefficient B defined at Proposition 3.4. 
Then we have: 

Proposition 4.5 

If conditions 1-6 and BwT,eq > BppAR-/-^eq 9.i"e satisfied, then 
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Biological prediction 8 (PPAR mutants) 

Under the hypotheses of Prop. 4.5 we have: 

1. PPAR knock-out reduces energy buffering. 

2. PUFA concentration increase under fasting is stronger in PPAR knocked-out cells compared to the 
same increase in wild type cells. 

Comment. The condition B^yr.eg > BppAfi_/_^^q is equivalent to Bi {R^f^)^^^^^+B2 {Rf°''^)^^^^^+ 
-^4^^ (iff) > 0, with Si > 0, -84^^ (iff) > 0. This means that even if B2 is negative 

the oxidation genetic control term is large enough to compensate. 

Experiments on transgenic mice showed that after a 72h-fast, fatty acids concentration increases at a 
higher extent in PPAR knocked-out cells with respect to wild type cells [BLC+04]. This is coherent with 
the observations by Lee et al.[La04] that for the same length of fasting time the hepatic accumulation 
of triacylglycerol is 2.8 fold higher in PPAR knocked-out than in wild-type mice. Hence, the global 
behavior of fatty acids is consistent with our predictions. 

Remark. Data from [La04] show a rather selective behavior among different hepatic TG PUFA in 
PPAR knocked-out mice: concentrations (in mg/g of liver) of a-Linolenic acid arc amplified on fasting 
to 4-fold higher levels than in WT mice, but Arachidonic, Docosahexaenoic and Eicosapentaenoic acids 
are depleted to non-detectable levels in PPAR mutants. Consequently, the behavior of some important 
regulating PUFA in mutants can not be explained by our model, at least not within the assumptions 
that we have made. In fact, our model is too crude to explain the contradictory behavior of part of the 
PUFA in mutants. A more complex model including more variables for PUFA and genes, could answer 
the questions raised by the experiment. Such a model should separate essential fatty acids from long- 
chain PUFA synthetized from the essential fatty acids, and include the genes involved in the synthesis 
of long-chain PUFA. 



4.4 Illustration of dynamics 

In order to illustrate the predictions of the previous section and the dynamical behavior of the model, 
we move here to another level of abstraction. The forms of the functions describing how primitive fluxes 
depend on concentrations are chosen, including numerical constants. The choice is rather generic, by no 
means precise. 

• Michaelis-Menten regulation functions. The characteristic features of those functions are the 
following: their value at is 0; their derivative is positive and strictly decreasing to at infinity. 
A typical analytical expression is: f{x) = -^f^, where k is the asymptotic value at infinity and 
k/K is the slope at the origin. This type of function models the dependency of an primitive flux 
on its substrate. 



• Repression functions. They are positive, strictly decreasing functions, tending to at infinity 
and have an inflexion point. A typical analytical expression is f{x) — j^^j^, where k is the 

tcla+l) ) 

a > 1 (also called Hill coefficient) is the logarithmic slope at the inflexion point. 



• Activation functions. They are positive, strictly increasing bounded functions and have an 
inflection point. A typical expression is /(a;) = b + j^^^a with K > 1. At 0, the value is b; the 

function tends to A; -|- 6 at infinity. K controls the position of the inflexion point: xe = ( ^^i^ ) 
a > 1 is the logarithmic slope at the inflexion point. 
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Figure 4.1: Starving/refeeding protocol, a) Concentrations of metabolites and main regulators. Starving 
begins at i = and refeeding at t = 750 (G = for < t < 750 and G = 10 for < < 0, t > 750). 
Fatty acids (Fi S/MU-FA, F2 PUFA) increase at fasting, with overshoots, decrease at refeeding, with 
undershoots. Acetyl-coA A and T (ATP, energy) have the opposite behavior. By PPAR and LXR we 
mean the active forms of nuclear receptors. These stimulate oxidation and ketonic exits at fasting (LXR 
falls, PPAR raises) and synthesis at refeeding (LXR raises, PPAR falls), b) Fluxes. At fasting (G = 0) 
there is practically no synthesis, and the Krebs cycle functions at a very low level. The fatty acid intake 
fluxes Fin are positive, oxidation and ketonic exit are strong. Refeeding inverses the Oxi,Kout/Syn, Krebs 
flux ratios, and changes the sign of the fatty acids intake fluxes (which become negative). 



• Energy dependence functions. Some reactions start when energy is missing and rapidly decrease 
to zero when energy increases. Other, like synthesis, are energy stimulated. We suppose that 
energetic regulation is sigmoidal, rather than hyperbolic. We use the following repressing function 
f{x) = j^^^i and the following activating function g{x) = . 

• Degradation functions. All products are supposed to degrade at a rate proportional to their 
concentration. 

• Consumption of ATP. Concerning ATP, we suppose that the consumption is linear. 

• Intake/ outake of fatty acids. The case of intake/outake of fatty acids is special. Both reactions 
are reversible. We suppose that outake is proportional to the internal concentration of fatty acids 
and intake occurs only when energy is missing inside the cell. 

Equations for metabolic fluxes 

Gly(G,T) = x^TfwTfe^ Krebs(A,T) = .^^^ Syn(A,T,EO = fe^I^rS- 

Oxi2(F2,T,E3) = ,:^^||^ Oxil(Fi,T,E2) = ^^i^|^^ DegT = ^tT 
Finl(Fi,T) = -fcpiniFi + ^^;^;"^T^ Fin2(F2,T) = -fcFin2F2 + x7|jTW Kout(A, E4) = ^^"^f^A 

Equations for genetic variables All controls are positive activations except the control of F2 on L: 
The final system in given in Table 4.1. 
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Figure 4.2: a) Response curves for varying food (glucose G). Quasi-stationarity (qs) corresponds to 
passing from a conventionally chosen normal, fixed feeding level G = 10 to a smaller, variable value of 
G, with enzyme levels unchanged from the normal feeding equilibrium. Genetic effect on energy balance 
: the energy T is lower at quasistationarity (no genetic regulation). Fatty acids increase with fasting 
(decreasing G) ; this behavior is more pronounced at quasistationarity. b) The strong lipolytic condition 

meaning that —gf- < has been plotted as a function of T, and for various values of G (0 < G < 10). 
c) Mode commutation; when food increases, the model commutes from a oxidation dominated mode to 
a synthesis dominated mode. 



Simulation of fasting/refeeding protocols In Fig. 4.1 we have simulated a fasting/refeeding protocol. 
One can notice the increase of fatty acids on fasting with an overshoot (as predicted in Prop. 4.2 and 
its Cor. 5). 

The dynamics has two timescales: a quick increase up to the maximum, then a slow decrease. The 
concentration of Acetyl-coA is decreasing at fasting. The energy (ATP) has an abrupt fall, then it 
recovers slowly as a result of oxidation. The behavior of nuclear receptors correspond to what is experi- 
mentally observed [Jum04]: LXR diminishes and PPAR is amplified at fasting. This leads to variations 
of the enzymes: oxidation E2,E3 and ketone exit enzymes E4 are amplified, the synthesis enzyme Ei is 
diminished. Again, this fits with experimental observations [Jum04]. Furthermore, fluxes have textbook 
behavior [Sal99] . Oxidation and ketone exist occurs during fasting, while synthesis occurs during normal 
feeding. During fasting fatty acids enter the cell Finl,Fin2 > 0. At normal feeding, Finl changes sign 
(de novo synthesized fatty acids exit) and Fin2 vanishes. The ketone overshoot at refeeding can be 
explained by enzyme inertia. The high fasting level of enzyme E4 can not drop immediately. As large 
amounts of Acetyl-coA are again available, this boosts the ketone production. 

The reader should be warned of a possibility not studied in this paper. Two timescales dynamics, 
involving the observability of the quasistationary states and the rapid overshoots and undershoots in 
the dynamics of metabolites and fluxes, can be avoided. For instance, glucose homeostasis could be 
responsible of slow instead of steep decrease of glucose at fasting. A slow input signal will drive the 
system quasi-statically, avoiding quasistationarity and rapid transients. As we do not study glucose 
dynamics and homeostasis, our model can not tell how glucose behave in time. This information should 
be provided by experiment. 

In Fig. 4. 2 we have simulated response curves of various metabolites, fluxes and energy when food G 
is changing. It can be noticed that fatty acids concentrations decrease with food and that Acetyl-coA 
concentration increases with food. Energy T is increasing with food. There is a buffering effect, that 
preserves energy against variations of food: energy T is not zero when food G is zero. As discussed in 
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Figure 4.3: Starving/refeeding protocol: dynamics of metabolites, energy and fluxes for a PPAR-/- 
mutant compared to wild-type. Fasting starts at i = 0, refeeding starts at t = 750. 



Prop. 4.4, a strong buffering effect means a weak slope of the dependence of T on G. Genetic regulation 
decreases this slope, therefore increases buffering as can be seen by comparing the curves in Fig. 4. 2a) 
at quasi-stationarity and at equilibrium. The antagonistic relation between synthesis and oxidation is 
illustrated in Fig. 4. 2c): when food G decreases, the synthesis dominated regime changes to an oxidation 
dominated regime. 

Satisfiability of the uniqueness and strong lipolytic conditions Our sufficient uniqueness condi- 
tions, allow, with small computational effort, to check the uniqueness of equilibrium for a given set of 
parameters. We have checked the strong lipolytic response condition for various sets of parameters. We 
have noticed that this condition is robust (see Fig. 4. 2b)). 

Simulation of the effect of PPAR mutations We have modeled a PPAR mutant by considering 
that the enzymes E2, E3, E4 controlled by PPAR have constant, unadjustable values. We have considered 
that the values of E2, E3, E4 are those for the normal feeding equilibrium state in wild type hepatocytes. 
fn Fig. 4. 3 the dynamics during a fasting/refeeding protocol are compared in the mutant and wild type 
case. The main feature of mutants is the difficulty to recover energy at fasting (see Fig. 4. 3 b)). This 
is the consequence of inefficient oxidation (notice the low oxidation fluxes in Fig. 4. 3 c)). Similarly, in 
mutants the ketone production is decreased ( see Fig. 4. 3 d)) and the fatty acids increase at fasting is 
more pronounced just like we have predicted qualitatively in Prop. 4.5 and its Cor.8(see Fig. 4. 3 a)). 



5 Discussion 

Let us first summarize some of the characteristics of our model. 

• ft is an integrative model, because it takes into account all the main processes of carbohydrate 
and fatty acid metabolism in liver (glycolysis, lipogenesis, Krebs cycle, fatty acids mobilization, 
oxidation, ketogenesis) together with their various regulation (metabolic, genetic). 

• Our model is not distributed: dynamical variables cope for average values in a tissue and no space 
information is taken into account. 

• The model that we propose is a low complexity abstraction. It has just enough complexity to 
represent the basic features of metabolism in the main nutritional states. Whenever possible, 
complex metabolic chains of reactions were modeled as a single global reaction preserving the 
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Table 4.1: A generic example for the model of the regulated metabolism. 



overall balance of products and reactants. Our main concern was to keep the model as qualitative 



• Our predictions are not dependent on specific numerical values of kinetic constants, on specific 
reaction mechanisms, or on specific forms of the functions relating fluxes to concentrations. We 
rather replace this information by sufficient qualitative conditions that are chosen as biologically 
significant as possible. 

• A qualitative approach has been used to discuss response properties. The introduction of the two 
types of states (quasi-stationary for rapid response and equilibrium for slow response) allowed us 
to distinguish between quick metabolic and slow genetic response. 

Our model copes with the main experimental findings on the behavior of regulated fatty acid 
metabolism in hepatocytes. Under fasting, the model shifts from a synthesis dominated regime to 
an oxidation/lipolysis dominated regime. This shift stabilizes energy, replacing food supply by reserve 
consumption. At short times, the shift is performed by metabolic control of synthesis, lipolysis and 
oxidation. At longer times, the regulatory effect of an increase of intracellular PUFA on the nuclear 
receptors PPAR and LXR reinforces this control. Refeeding shifts the system in the opposite direction. 

Our model is sufficiently general to apply to various higher organisms. Nevertheless, there are some 
biases in this model. More precisely, we have considered only hepatocytes, which is justified by the fact 
that in various species like chicken, rodents and humen, fatty acid synthesis and oxidation occur in the 
same organ (liver). This is not so for other species such as pigs, for which synthesis takes place mainly 
in adipocytes. 

We have also proposed a methodology to build small complexity abstractions that integrate various 
qualitative aspects of regulated metabolism. These abstractions are by no means rigid. On the contrary, 
they are evolutive and can integrate new experimental results. As an example of possible evolution of 
the model presented here we should mention the role of various PUFA, already discussed. 



From here on, we develop the mathematical setting and arguments which lead to the results announced 



as possible. 



6 Method 



before. 
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6.1 Switches, shifts and obervabiUty of equilibria 



Equilibrium switches and shifts. The two types of equilibrium changes shift and switch can be 
mathematically described as follows. Let Z = (X,Y) and F(Z,p), Z e [/ c M", p € J C M, be a 
differentiable vector field defined on an open set U of M" depending smoothly on a parameter evolving in 
an open interval /. We suppose that for any p in a subinterval J C /, the vector field i^p : Z ^ F{7i, p) 
admits at least a singular (equilibrium) point, that is: there exist Z € U such that Fp{Z) = 0. We call 
the 0-level (equilibrium) curve: 



£ is a differentiable curve in J x ?7. According to its shape, we get a switch or a shift. More precisely, 

1. A shift is characterized by the following features: 

• £ is a graph Z = L{p). Equivalently it means that for each value of p G J there exist a unique 
equilibrium point, 

• £ is a sigmoid. That is there exist a threshold value po, an interval ]po — 5, Po + S[, two 
equilibrium values Zi and Z2 and e > such that for p < po — S, \L{p) — Zi| < e and for 
p>Po + (5, |L(p)-Z2| <e. 

See an illustration in Fig. 6.1. 

The sigmoid shape of the level curve has two consequences: a "jump" effect and reversibility. If 
we start with a certain value of the parameter p/, say p/ < po — S, and increase smoothly the 
parameter up to a final value p^? > po + (J, we first observe an asymptotic state close to Zi and a 
sudden "jump" to an asymptotic state close to Z2. We insist on that the apparent "jump" is due 
to the steepness of C at po which induces a small 6: in reality there is no discontinuity. Now if we 
start with the value pp and decrease smoothly the parameter, we note exactly the reverse: a jump 
at po from an asymptotic state close to Z2 to an asymptotic state close to Zi, a property so-called 
reversibility. 

Notice that both properties are often used to identify shift like phenomena. 

2. A switch has the following properties: 

• £ is not a graph. In particular, there exist two parameter values po and pi for each of which 
the vector field has two singularities (equilibria) and for any p e]po,Pi[, the vector field Fp 
has three singularities (one unstable and two stable). For p < Po or p > pi there is only one 

singularity (stable). 

• Po and pi are bifurcation points of saddle node type, that is an attracting and a repelling 
singularity collide and disappear. 

• There exist S > 0, two equilibrium values Zi and Z2 and e > such that for p < po — 5, 

[L(p) - Zi| < e and for p > pi + (5, \L{p) - Z2I < e. 

To understand the hysteresis effect implied by the shape of the curve £ it is worth considering 
the experimental curve, that is, the curve of observed equilibria when moving the parameter. If 
we start with a certain value of the parameter p/, say p/ < po — 5, and increase smoothly the 
parameter up to a final value Pi? > po + S, we first observe an asymptotic state close to Zi and a 
sudden "jump" to an asymptotic state close to Z2. To the contrary to the previous shift situation, 
this jump is a real discontinuity. Notice that the jump occurs at the parameter value pi. Now 
if we start with the value pp and decrease smoothly the parameter, we note a jump at po from 
an asymptotic state close to Z2 to an asymptotic state close to Zi. We do not have reversibility, 
because the jump occurs for different critical values of p. The curves in both cases are depicted in 
Fig. 6.1. 

Experimental curves and equilibria. Consider now the differential system: 



£ = {(p,Z) eJxU,Fp{Z)=0}. 



(6.1) 



— = F(Z,p), ZeRl,pelcR 



(6.2) 
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Figure 6.1: Switches and shifts 



The points on the 0-level curve of the vector field F are equilibria of the dynamical system (Eq.(6.2)). 
In the following we consider the simple situation when for any p G I there is a unique attractive equi- 
librium Zo(p), of open attractive basin B{p). Let us recall that the attraction basin consists of all Z 
such that the trajectory u{t,Z,p) of the system (6.2), starting at Z, approaches Zo(p) for large times 
\u{t,Z,p) -Za{p)\ ^ (t-^oo). 

A typical experimentation consists in starting in an equilibrium Zo(pi) and suddenly changing the 
value of the parameter p from pi to p2- Suppose that the following condition is fulfilled Zo(pi) G B{p2). 
This condition is automatically fulfilled if for instance Zo(p2) is globally attractive, i.e. B{p2) = M". 
Then, the observed state u(t, Zo(pi),P2) will approach the state Zo(p2) on the 0-level curve after a long 
enough time. Suppose now that the 0-level curve Z = L{p) is such that the component Li{p) satisfies 

> for any p £ I. Then, it exists T such that Ui{t, Zo(pi),P2) > (^o)j(pi) for t > T, meaning that 
we observe an increase of the component i of the state of the system between the beginning and the end 
of the experimentation. 

Hence, the shape of the 0-level curve informs on variations of products during an experiment. This 
justifies that we can predict experimental behaviors of the system from the study of its equilibria. 

Observability of quasi-stationarity. Let us suppose now that the system (6.2) has two time scales, 
one slow and one fast. This can be taken into account by supposing that there is a small parameter e 
representing the fast time scales and that the system (6.2) reads ^ = $(X, Y,p), ^ — e'4'(X, Y,p), 
where X £ R"^, Y £ R"% Uf + ns — n. Considering the time scale t ~ et we arrive to the more classical 
form: 



A result due to Tikhonov [TVS80] and reformulated by Fenichel [Fen79] implies that the system (6.3) 
has a remarkable behavior. Suppose that the following two conditions are satisfied: 

• For fixed Y = Yq, the system ^ = *(X,Yo,p) has an attractor Xo{Yq,p) that satisfies the 
relation '>I'(Xo(Yo,p), Yo,p) — 0. The equation X = Xo(Yo,p) defined the slow manifold Sod?). 

• The Jacobian matrix _Dx admits n j eigenvalues with strictly negative real parts at (Xq ( Yq , p) , Yq ) . 

Then the system (6.3) admits an attractive invariant manifold I](e,p) close to the slow manifold T,o{p). 
Thus, trajectories sufficiently close to the slow manifold converge to it. The smaller e is, the quicker is the 
convergence. Notice that the slow manifold defines the quasi-stationary states. The Tikhonov-Fenichel 
result means that the system rapidly tends to quasi-stationarity. 

6.2 Existence of equilibrium (Theorem 3.1) 

The proof of Theorem 3.1 uses the following more general theorem: 
Theorem 6.1 

Let $(X) = G(X) — A(X) be a smooth vector field on R" (M" represents all the vectors of M" having 
non-negative coordinates) such that : 




(6.3) 
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1. G is bounded, 

2. For aJJ X = (Xi, . . . , X„) such that Xi = and Xj ^ for all j ^ i, G satisfies Gi(X) > 0, 

3. A = (Ai(Xi), . . . , A„(X„)) : Wl — *■ W^, and Aj are differentiable and satisfy Aj(0) = and 
lim||x||->+cx3 Ai(X) = +0O, for all 1 < i < n. 

Then the equation *(X) = has at least one solution in R" . 

The proof of the Theorem is based on the following standard mathematical lemma. 
Lemma 6.1 

Let D be a smooth ball in R" and let S be the boundary of D. Let $ be a differentiable vector field 
defined on a neighborhood of D. If # points inward D at any point of S then * admits a zero in the 
interior of D. 

Sketch of the proof of Lemma 6.1. By the Poincare-Hopf formula a sufBcient condition for having 
a zero in the interior of D is to have a non-zero index for the vector field on S. Since $ points inward 
D on S, we can construct a smooth change of variables which conjugates * on a neighborhood of D to 
a vector field defined on a neighborhood of the unit n— ball B", such that on a neighborhood of the 
unit n— sphere S", coincides with the radial vector field X —X. For this vector field we can 
compute its index, which is 1 or —1 according to the parity of n. The Lemma is proved since the index 
is a differential invariant. ■ 

Proof of Theorem 6.1. Prom Lemma 6.1, it is enough to find a smooth ball in the positive orthant 

on the boundary of which the vector field $ points inwards. 

For R > 0, let us consider the intersection domain of the closed n-ball of radius R with the positive 
orthant: A = {X G M" , ||X|| < R}. This domain is a topological ball; let us denote S its boundary. If 
X e E and none of its components is 0, then for R large enough, *&(X) points inward A, because G is 
bounded and Aj(X) tend to infinity with X, hence $i(X) < 0, for all 1 < i < n. On the other hand, if 
only one of the components of X is 0, then by hypothesis (2), #(X) points inward A. Since the set of 
points where the property of pointing inwards is open, we can find a smooth ball D contained in A and 
sufficiently close to it, such that on the boundary of D, the $ points inward D. ■ 

Actually, the Theorem implies naturally a stronger result which will be useful in the proof of unique- 
ness of equilibrium. 

Corollary 6.1 

Under the hypotheses of Theorem 6.1, Jet X = (Xi,X2) be any partition of the variables. We write 
<l>i(Xi, X2) for the projection of <I>(Xi, X2) in the vector space spanned by the coordinates of Xi. Given 
X2, the system of equations <l?i(Xi,X2) = 0, where X2 is considered as a constant parameter vector, 
admits a solution in Xi with non negative entries. 

Proof. Theorem 6.1 applies to the vector field $i(Xi,X2), where X2 is considered as a constant 
parameter vector, since $1 satisfies the hypotheses of the Theorem 6.1 as soon as $ satisfies them. ■ 

Proof of Theorem 3.1. The proof runs by applying Theorem 6.1 to the following vectors: 

/ SyUpgq - Oxilpeq + Fiul \ 

-Oxi2 + Fin2 
acGly — asSyn + axKrebs + aoiOxil + ao20xi2 
y Gly + niOxil + 7x2 0x12 — Krebs — Kout — niSyn J 

Notice that equilibrium states of the implicit models are the zeroes of the vector field ^'(X) = 
G(X)-A(X). 

Let us verify hypotheses of Theorem 6.1. First, G is differentiable because all fiuxes are differentiable. 
Then, G is bounded because it is composed of primitive fluxes which saturate at high concentrations of 
metabolites (Condition 3). Finally, to verify the condition Gi{. . . , Xi = 0, . . .) > 0. it is enough to notice 
that each coordinate can be decomposed into the difference of the fltixes which produce the variable 
and the fluxes which consume the variable. The second ones are zero when the variable is zero by zero 
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substrate effect (Condition 3). The sum of the producing fluxes is strictly positive by recovery effect 
(Condition 3). 

The assumptions on A are satisfied by Hnearity of degradation terms and by unboundedness of ATP 
consumption term (Condition 3). ■ 

6.3 Box reduction of systems of non-linear equations 

Let $ : R" x A ^ R", where A is a compact subset of M"^, be a differentiable vector field. $ defines the 
following system of linear equations parametrized by p: 

>S:*(X,p) = (6.4) 

Box of a system of equations. We call box of the system (6.4) a subset X^*) of the set of variables X, 
such that X = (X^*) , X''^' ) is a partition of the set of variables. The variables X^*) , X^^^ are called internal 
and external variables, respectively. To each partition of the variables, let us consider the corresponding 
partition of the vector field components # = (#(*\ $(*)). 

We call box equilibration the elimination of internal variables from the equations defined by the 
internal part of the vector field: *(*)(X(*), X(«),p) = 0. 

Sequence of box equilibration. After a box equilibration the internal variables can be expressed as 
functions of the external variables. A sequence of box equilibrations is the finite iteration of the following 
operations: 

1. Define Xi = X and *i(Xi,p) = *(X,p), X>i(p) = R". 

2. At /c-th iteration, divide the variables and the vector field components into internal and external 
parts X, = (X«,X(^)), *. = (*«, *(^)). 

3. If the external part is not empty then: 

• solve <&^*''(X^*\ X^"^-*, p) = with the constraint (X^*\x^*^^) e I'a;(p) and express the internal 

variables as functions of the external variables X^'^ = M.]^{X.^\-pi). Notice that the solution 
might not be unique, that is A^k is not necessarily univocal. We restrict our discussion to 
the case when the number of solutions is finite and bounded, such as for polynomial systems. 
Also, notice that one has a solution for X^'^'' in a maximal domain Vk+i{V')- If 2^fc+i(p) is 
empty then stop: there is no solution. 

• define Xk+i = X^f\ and *k+i = *^r'(A1k(xL^\p),xL^\p). 

4. If the external part is empty then solve *[,'^(Xj.'\p) = and stop. Conventionally, in this case 
■Dfe+i (p) is considered non-empty iff the equation has a solution. 

5. go to step 2. 

A sequence of box equilibrations is complete if all components are equilibrated i.e. 

x = xf ®X^'^ ©...©X^^^. 

After a complete sequence of box equilibrations one should be able to express metabolite levels as 
functions of the external parameters: X = Al(p), where M. results from a composition of the functions 

An example of (incomplete) sequence of box equilibrations is the reduction (at equilibrium) of genetic 
variables in the mixed metabolic/genetic differential system. There is only one box whose internal 
variables are the genetic variables Y. These are eliminated from the equations -(/"(XjY, p) = 0. The 
functions A^i(X,p) are YP'=9(X.p), and the reduced fluxes are *2(X,p) = *(X, Yp«9(X, p), p). 

Existence and uniqueness of solutions. 

Box equilibrations perform nothing else than the substitution method for non-linear systems of equa- 
tions. The existence and properties of solutions relatively to box equilibrations are straightforward. 
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Proposition 6.1 

• A solution of the system (6.4) exists for a value of the parameter p if there is a complete sequence 
of box equilibrations with non-empty domains T>k+i{p). 

• The function M is univocal (to one p corresponds a single value of M.) if all the domains Vk+i (p) 
arc non-empty and each one of the function Alk is univocal on its maximal domain Dfc_|_i(p) for a 
complete sequence of box equilibrations. 

This property is useful to prove the existence and uniqueness of solutions of systems of non-linear 
equations. It is enough to choose a complete sequence of box equilibrations and to show that at each 
step the functions M.k are univocal on non-empty domains Vk+i{p). 

It is difficult to give a "only if version of the property. Indeed, even if we find a box such that the 
equations (X^*\ x[,'^\ p) = have multiple solutions in X^'^ it is not excluded that some of these 
solutions are incompatible with the rest of the equations: after all the box equilibrations we may still 
have an unique solution. 

Sketch of the proof of Proposition 3.2 and Theorem 3.2. According to Proposition 6.1 a sufficient 
condition for uniqueness of equilibrium is to find a complete sequence of box equilibrations for the 
state equations of metabolic variables (Eq.(2.3)). To simplify notations, in this section we write the 
reduced flxrxes without the subscripts peq, gnr. Notice that in order to get the state equations a ffist box 
elimination of the genetic variables (shown to be univocal) has already been performed at equilibrium. For 
genetically non-regulated states, the reduction of the genetic variables has been performed by replacing 
them by constants. 

Let us consider a box decomposition corresponding to the following partition of variables: 

• Box 1: {A,Fi,F2}. 

• Box 2: {T}. 

The main steps of the proof are the following: 

• Step 1: the equilibration equations for Box 1 with internal variables A, Fi, F2 is: 

$a(G,A,Fi,F2,T) = (6.5) 
$Fi(A,Fi,F2,T) = (6.6) 
$F.(F2,T) = (6.7) 

In Lemma 6.3, we prove that for all (G,T) G [0, Gmao;] x this system has a unique solution 
(A(i)(G,T),Fi(i)(G,T),F2^^)(G,T)) e ^ which ends the proof of Proposition 3.2. 

• Step 2: the equilibration equation for Box 2 with internal variable T is: 

$t(G,A(i)(G,T),FiW(G,T),F2^^^(G,T),T) = (6.8) 

In Lemma 6.4, we shall provide a sufficient condition for this equation to have a unique posi- 
tive solution T = T(^)(G'). The unique solution of the state equations is thus A = PS'^\G) = 
A(i)(G,T(2)(G)), Fi = Fi(2)(G) = Fi(i)(G,T(2)(G)), F2 = F2(')(G) = Y2^^\G,1:'^^^{G)). This 
solution depends on being at equilibrium or at quasi-stationarity, because the expressions of the 
reduced fluxes $a, ^Fi, depend on being at equilibrium or at quasi-stationarity. 

• Step 3: finally, we show that the strong lipolytic condition is the sufficient condition of Lemma 6.4. 

All the formal manipulations of this sequence where performed using Wolfram Research Mathematica 
version 5.2 software. 

To treat the case of Box 1, we use the following result which is a direct consequence of Gale-Nikaido- 
Inada theorem [Par83] . This theorem can be seen as a generalization to higher dimensions of the mono- 
tonicity of functions on R. Let us recall that a principal minor of a matrix M = (?Tii,j)i,je{i,...,n} is 
defined as detM/, where / C {1, . . . ,n} and Mj = {mij)ij^i. 
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Lemma 6.2 (Gale-Nikaido) 

If {x,y,z) {fx,fv,fz) is a diffcrcntiable mapping from M.'^ to M.^, of Jacobian J, such that all the 
principal minors of — J are positive, then this mapping is globally univalent. In particular the system 
fx{x, y, z) = 0, fy{x, y, z) = 0, fz{x, y,z) = has a unique solution if a solution exists. 

Lemma 6.3 

For all (G, T) e [0, Gmoa;] x R+, the system of Eqs. (6. 5), (6. 6) and (6.7) has a unique solution in 
(A, Fi,F2) e M.^, both at genetic partial equilibrium and at fixed genetic variables . This expresses 
A,Fi,F2 as univocal functions of G and T, denoted by A^-^^CT), Fi(^)(G,T) and F2^^^(G,T). 

Proof. We apply Corollary 6.1 to prove that the system of equations (6.5), (6.6) and (6.7) has a solution 

for every fixed (G, T). 

Let us consider the mapping (A, Fi,F2) ($a,^Fi,$F2) on M.^. To prove the uniqueness of the 
solution to the system, we apply Lemma 6.2 to this mapping by ensuring that all the principal minors 
of — J(^^ are positive, where J^^^ is the Jacobian of this mapping: 



/_V*°* ^'I'A \ / tot „ ^,Oxil „ pOxi2 , „ nOxil , „ pSyn nKoul 

AA aFi dF2 \ (-Xa ^IXFi '^2-Kf2 +"l-n-F2 + '^1^F2 ^ ^F2 

'9'i'Fi ,^iot Q^Fi ^ Syn tot pOxil pSyn 

~aA~ AFi "OTT I AA AFi -"■F2 -"•F2 

V -x^°7 V -x'fI 



The principals minors of -J^^) are aU positive: Xa* > 0, Xa*Xf* - "iXA^'Xpr'^ = Xa*Xf*(1 - 
^Syn^Oxii) > 0, Xf°*Xa*Xf°*(1 - PA^'^P^f^) > 0, as a consequence of Prop. 3.3. This is valid both at 
genetic partial equihbrium and at fixed genetic variables . 

Hence, Lemma 6.2 applies, so that there exist functions A(i)(G,T), Fi(^)(G,T) and F2^^^(G,T) that 
arc the unique solutions of the system of Eqs. (6.5), (6.6) and (6.7) for each (G,T). These functions arc 
difFcrcntiablc on by the implicit function theorem. ■ 

Proof of Proposition 3.2. Proposition 3.2 derives directly from Lemma 6.3. ■ 
Lemma 6.4 

Eq.(6.8) has a unique solution in T as soon as the following inequality holds for every (G, T) g [0, Gmax] x 
M+: 

^+(aKXr^^-a5Xr)^+aoixFr^^ + («oii??2^'^+ao2iiFr+asi?r)^ < 0- (6.9) 

Proof. Let $t'^^(G,T) = $t(G, A(i)(G,T),FiW(G,T),F2(^)(G,T),T). The biological hypotheses 
imply that Theorem 3.1 applies that is, for every G, there exists (a, /i, /2, t) G such that 

$A(a,/i,/2,i) = $Fi(a,/i,/2,t) = $F2(/2,i) = $T(G,a,/i,/2,t) =0. 

By uniqueness in the previous Lemmas, we get a = A(^)(G,t), /i — Fi'-^-'(G,t) and /2 — F2^^^(G,i). 
Hence ^t'^^G,*) = and the function $T*^^(G,r) has a root in T for every G. 

Moreover, the function $t^^^ is differentiable on R^. From the definition of the function it 
follows: 

—q^ = -qY + {o^kXa -asXA )-^+aoiXFi grj. +(qoiJ^F2 +ao2i?F2 +Qs-Hf; ) (6.10) 

Hence if the inequality 6.9 is satisfied, then '^^'!p is negative. In other words, ^x'^' is monotonic so 
that it has a unique zero. ■ 

Proof of Theorem 3.2 Theorem 3.2 derives directly from Lemma 6.4. ■ 
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Lemma 6.5 

The derivatives , ^gij , ^^p— can be expressed by means of fluxes and of control coefRcients in 
the following way: 

det(J^^)) aA'^) 



tot f ryfjrly T->Krebs . _ r Oxi2 7->Fin2 . / Oxi2 \ 7-,Oxi2-] , r/-, Oxil \ / j-jOxil . 7-,Syn\ . 

Xf2 I-^t - -^T +n2|pF2 + (1 - Pf2 )Rt i + - Pfi )(^t -\-Rt ) + 



tot 



Xf7 





Oxil oFinlii . / 7->Fin2 7->Oxi2\r oKout , „ / oO^il i DSyn\/-| Oxilxi 

pFi -Rt J} + (-Rt - Rt I[~Rf2 +ni{RF^ +i?p;)(l-pFi )J 
_det{Jp_m^ = X^-.Hpr - ^i?'^"^ + n^Rr') + ni [ii?'"^ - (1 - + ^1"")]} 

+{Rt - Rt )[ni{RF2 + -Rps ) - «2-Rf2 + p/ Rf2 \ 
£^ = (6.11) 

wiere - det(j(i)) = x'fIxTx'fH^ " pI'VR"'')- 
Proof. The lemma follows straightforwardly from 



Lemma 6.6 

= [A{R^-'' - + BiR^-'^ - + c- D]/[m{l - pfy^f')] (6.12) 

where A, B, C, D are combinations of control parameters defined in Eq. (3.2). 

Furthermore, if the stoechiometric condition 5 is fulfilled, then X > 0,^4 > 0,-Bi > 0,-64 > 0,C > 
0,0 > 0. 

Proof. The proof is a lengthy but straightforward formal manipulation of Eqs.(6.10),(6.11). We have 
gathered control coefficients into as large as possible positive combinations. 

As an illustration of how the stoechiometric condition was used let us consider the sign of D. From 
aoi > as, X/m - aoi(l - pi^") = (aoi - as)pi^" + maKpA'''''" > 0, p^/" < 1, it follows that D > 0. 

■ 

Proof of Proposition 3.4. It follows directly from Lemma 6.6. ■ 
Lemma 6.7 

Let T = T(^'(G) be the solution of the equation (6.8). If the strong lipolytic response condition and the 
stoechimetric conditions are satisfied then '^'^^ > 0. 

Proof. The partial derivatives with respect to G of the metabolites after equilibration of the box 1 can 
be obtained from: 






tot \ 



TjGly / XFi 
, Syn 



I ^tot^tot(-i _ „ayn Oxil\ 

/ Xa Xfi Pa Pfi ) 



From this and from the definition of it follows 

(6.13) 

If as < aoi < niac (which is true by the stoechiometry condition 5), then < aoi(l — Pa'") = 

r^,^ ^^ ^Oxil^Synx , .Oxil _Syn Syn /lOxil Syns, , „Oxil .Syn ^„„Syn , 

aoili - Pfi Pa ) + Q^oiPpi Pa ~ Q^oiPa < ?^iQ:g(1 — Pfi Pa J + ctoiPpi Pa ~ oisP\ , ana 
consequently ' > 0. 
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The strong lipolytic response condition is equivalent to having ' < 0. If this is satisfied, then 

dG — \ dG J ' \ dT J ^ ^- ■ 

Proof of Proposition 4.3. It follows directly from Lemmas 6.7. ■ 
Lemma 6.8 

Let F^^^(G) = F^^^(G,T(2)(G)). If the strong lipolytic condition is satisfied, then the sign of is 

equal to the sign of R^""'^ - 

Proof. The chain rule gives 

dG ~ dG ^ dT dG ' 

From the proof of Lemma 6.7, ' = 0. It follows from Lemma 6.5 that the sign of ^ ' is the 
same as the sign of /i^'^'^ — Moreover, if the strong lipolytic condition and the stoechimetric 

condition are satisfied, then > 0. ■ 

Proof of Proposition 4.1. It follows directly from Lemmas 6.7,6.8. ■ 

Proof of Propositions 4.2 and 4.4. The differences between equilibrium and quasi-stationarity occur 
at two levels: 

1. At quasi-stationarity F2 does not regulate the genetic variables: 

2. At quasi-stationarity the control of F2 on its oxidation is only a metabolic substrate effect. Genetic 
control is added at equilibrium. We have R^^^ = + fff with > 0, and Xf* = 

^Oxi2 _ OF^^ ^o.i2 ^ (1 _ ^ith ^ < 0. Furthermore, (fff )^ > and (f||) 



eq 

Hence: 

)e, > > ix^De, > {x'fD,, , {P9f)^, > {P9f),, ■ (6.15) 

We deduce from Eq.(6.13) that ^"^^^ ' is the same at equilibrium and at quasi-stationarity. Prom 
Lemma 6.6, it follows : 



-^)^ - ^ ~ „i(i-pf-p°-i) - -Rt""), (6.16) 

where 7?. is a term not changing from quasi-stationarity to equilibrium and the expression of B is detailed 
in Prop. 3.4. 

Prom = -(^^^) /{^^^), from Bg, < B^g, and from Eqs.(6.14,6.15,6.16), it follows 

that f^^) > (^) . From ^ = -^(i?Fin2 _ B°^^)/^ta ^nd Eq.(6.15) it follows that 

\ /as V / eo 



.(2) 



dG 



> 

qs 



.(2) 



dG 



eq 

Proof of Proposition 4.5 We follows closely the proof of Prop. 4.2. The differences between PPAR - 
/— and WT cells occur at two levels: 



(ijo-i) = f i?Fr) = (6.17) 

V ^ J PPAR-/- V 2 I ppAp_i_ 



PPAR-/- \ ' y PPAR-/ 
( r)Oxi2\ - / rjOxi2\ / tot\ . I tot\ 

y^^' )wT,e, > JPPAR-/-,.,' (^f2 VT,e, > (Xf J pp^«_/_,,, , 

) > (p'if ) (6.18) 

V ^ /WT,eq \ ^ J PPAR-/-,eq 

If BwT,eq > Bppj^ji_/_ gg, it follows (along the same lines as the proof of Prop. 4.2) that 
{'^)pPAR-/-,.q > ('^)^T,e.- ^ = -'^(^"" - ^"')/XF* and Eq.(6.18) it follows that 

I '^'^ \eq,PPAR-/- I '"^ \eq,WT' 
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